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ABSTRACT 


A  user  friendly  program  was  developed  to  calculate  fadgue  life  using  Strain- 
Life  equations,  given  either  a  stress  or  strain  history.  Additionally,  the  material 
parameters  and  associated  stress  concentration  factors  can  be  varied.  Since  certain 
material  constants,  such  as  cyclic  strength  coefficients  (K’)  and  strain  hardening 
exponents  (n’)  vary  during  a  material’s  fatigue  life,  the  program  is  capable  of  either 
keeping  them  constant  or  varying  them  as  a  function  of  elapsed  cycles.  The  program 
was  then  utilized  to  examine  the  effects  of  varying  K’  and  n’  on  the  calculated  fatigue 
life  of  aluminum  7075-T6  under  a  typical  flight  load  history. 
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I. 


INTRODUCTION 


There  are  distinctly  different  approaches  used  to 
calculate  cumulative  fatigue  damage  during  the  crack 
initiation  and  crack  propagation  stages.  Since  the  United 
States  Navy,  Naval  Air  Systems  Command  considers  any  cracked 
part  to  be  "failed",  the  focus  of  this  thesis  was  on  crack 
initiation. 

The  definition  of  fatigue  damage  during  the  initiation 
phase  of  a  crack's  life  is  difficult.  The  damage  during  the 
crack  initiation  phase  can  be  related  to  dislocations,  slip 
bands,  microcracks,  etc,  but  the  phenomena  are  microscopic  and 
are  not  easily  correlated  with  macroscopic  measurements. 
Because  of  this,  the  damage  summing  methods  typically  used  to 
calculate  crack  initiation  are  empirical  in  nature.  They 
relate  damage  to  life  consumed.  Life  here  refers  to  the 
physical  separation  of  a  small  test  specimen  which  is 
subsequently  used  to  approximate  crack  initiation  in  larger 
aircraft  components,  since  larger  components  will  have  a  much 
larger  critical  crack  length. 

The  most  common  method  for  summing  damage  is  to  use  the 
linear  damage  rule,  also  known  as  Miner's  rule: 


(1) 
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where  ni  is  the  number  of  cycles  applied  at  a  given  stress 
level  and  is  the  total  fatigue  life  for  constant  amplitude 
loading,  at  that  stress  level. 

When  applying  Miner's  rule  to  a  variable  load  history,  the 
life  used  up  at  each  load  change  has  to  be  calculated.  This 
can  be  done  with  one  of  several  possible  strain- life 
equations;  for  instance: 

Morrow's:  —  =  —  (2Nf)  b  +  €f  (2Nf)  c  (2) 

2  E 


Manson- Hal  ford: 
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(2Nf)b+e'A 


_C 

b 


(2  Nf)c 
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A  (  _/  \  2  (4) 

Smith- Watson-Topper :  (2N£)  2b+oWf  (2Nf)  b*c 

2  E 

All  three  of  these  equations  require  the  local  change  in 
strain  and  either  the  maximum  or  mean  local  stress.  These  can 
be  obtained  from  relating  the  far  field  stress  to  the  local 
stress  using  Neuber's  empirical  rule: 

KtSe  =  oe  (5) 

Relating  stress  changes  to  strain  changes  can  be  done  with  the 
monotonic  stress- strain  equation: 
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(6) 
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or  the  hysteresis  equation: 

1 

_Ae  _  da  _  Act  lTi7  (7) 

2  2E  2Kf 

[Reference  1] 

Programs  are  available  to  solve  strain- life  equations; 
however,  a  fatigue  life  program  was  needed  for  research  at  the 
Naval  Postgraduate  School,  where  the  user  can  modify  the 
solution  algorithms  or  have  access  to  the  program's  source 
code  to  make  changes  and  explore  various  facets  of  the  theory. 

For  example,  to  evaluate  the  effects  of  varying  the  cyclic 
strengch  coefficient  and  the  cyclic  strain  hardening  exponent 
throughout  the  duration  of  applied  loading,  a  program  named 
FLP  was  developed.  This  program  solved  the  three  strain- life 
equations,  (Equations  (2)  through  (4),)  using  either  fixed  or 
varying  values  for  n'  and  K'  as  specified  by  the  user. 

To  develop  a  realistic  sequence  of  loads  from  a  known 
spectrum,  LOADGEN  was  created.  It  was  used  to  develop  a 
realistic  load  history  for  an  A-6  based  on  the  three  a  "g" 
count  data  developed  by  LT  Rich  Walters  [Reference  2] .  The 
load  sequences  were  processed  with  the  program  using  various 
strain-life  equations. 


II.  STRAIN-LIFE  COMPUTATION  PROGRAM 


A.  GENERAL 

The  Fatigue  Life  Program  ( FLP )  was  written  using 
Microsoft's  QuickBasic.  QuickBasic  is  a  relatively  simple 
programming  language,  but  an  updated  and  more  capable  version 
of  BASIC,  which  allows  the  program  to  be  compiled  into  an  .EXE 
file,  executable  on  any  computer  using  MS-DOS  or  PC-DOS 
operating  system  2.1  or  later. 

FLP  was  designed  to  be  as  "user  friendly"  as  possible, 
and  to  the  greatest  extent  practical,  menus  were  used  to  set 
the  various  options.  The  program  was  built  in  a  modular 
fashion  and  documentation  was  included  throughout  the  program 
to  facilitate  debugging  and  later  modifications. 

B .  MAIN  PROGRAM 

The  main  program  controls  the  general  processing  flow.  It 
starts  by  establishing  certain  constants  and  variables  to  be 
used  throughout  the  program,  and  gives  initial  values  to  most 
of  the  user  selected  options,  which  could  be  considered 
default  values.  Examples  of  default  values  would  be  British 
units  and  aluminum  7075,  which  appear  as  being  the  units  and 
material  selected  when  the  program  is  first  started. 

The  program  then  enters  a  perpetual  DO  loop  wherein  it 
sets  all  the  user  definable  options,  reads  in  the  specified 


load  sequence,  processes  the  load  sequence,  and  outputs 
processed  data  to  the  designated  file  before  returning  to  the 
user  definable  options.  The  loop  is  exited  and  the  program 
terminated  by  pressing  ESCAPE  when  the  option  menus  are  being 
displayed. 

The  main  program  also  sets  the  video  configuration  for  the 
program,  allowing  it  to  automatically  select  the  best  mode 
available,  and  incorporating  error  trapping  sequences  to 
indicate  when  insufficient  memory  is  available. 

C.  USER  SELECTABLE  OPTIONS 

The  selection  of  the  various  options  available  to  the  user 
is  accomplished  through  several  menus,  each  of  which  is  a 
subroutine,  or  set  of  subroutines.  There  are  two  menus  to 
define  user  selectable  options,  and  a  third  to  input 
properties  for  new  materials.  All  three  of  these  menus  print 
a  list  of  key  instructions  at  the  top  of  the  screen,  and  call 
on  another  subroutine  to  update  the  display  screen. 

The  first  menu  to  appear  is  shown  in  Figure  1  and  allows 
the  operator  to: 

•  choose  between  using  British/U.S.  (Brit)  units  or  Standard 
International  (SI)  units; 

•  choose  a  local  stress  concentration  factor  (Kt)  ; 

•  choose  a  method  to  calculate  the  fatigue  stress 
concentration  factor,  either  by  use  of  Neuber's  method, 
Peterson's  method,  or  manual  entry; 

•  choose  another  screen  mode  if  more  than  one  is  available; 
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•  choose  a  material  from  the  existing  material  data  base  or 
manually  enter  a  new  material,  which  may  be  saved  in  the 
material  data  base; 

•  review  the  material  properties  associated  with  the 
selected  material. 

Prior  to  exiting  this  menu,  the  program  checks  to  ensure  that 
a  value  has  been  either  entered  or  calculated  for  the  fatigue 
stress  concentration  factor  (Kf) .  If  a  value  for  Kf  other 
than  zero  isn't  present,  the  program  will  prompt  the  user  for 
one.  The  program  may  also  be  terminated  normally  from  the 
first  options  menu  by  pressing  ESCAPE. 


UP  .  Move  to  next  field 

DOWN  .  Move  to  previous  field 

LEFT/RIGHT  ...  Change  field  up/down 
El  ....  Display  matrial ’s  parameters 
ENTER  ....  Start  with  current  values 
ESCAPE  .  Quit  Program 


Type  of  units  (S[  or  British)  [  Brit  ] 

Material  [  Aluminum  7075-T6  as  recieved  ] 

Stress  concentration  factor  (Kt )  [  1.000  ] 

Method  to  calculate  Kf  [  Manually  Entered  ] 

Fatigue  stress  cone,  factor  (Kf)  [  0.000  ] 

Screen  Mode  [  12  ] 


Figure  1  Options  Menu  l 
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The  second  menu  of  user  selectable  functions,  shown  in 
Figure  2,  appears  when  the  operator  exits  the  first  menu  by 
pressing  RETURN.  It  allows  the  operator  to: 

•  choose  between  using  stress  or  strains  as  input  loads: 

•  choose  which  strain- life  equation  will  be  used,  either 
Morrow's,  Manson-Half ord,  or  Smith-Watson-Topper ; 

•  choose  between  using  either  fixed  values  for  the  cyclic 
strength  coefficient  and  cyclic  strain  hardening  exponent 
or  using  values  that  are  a  function  of  elapsed  cycles; 

•  choose  between  calculating  the  cycles  to  failure,  the 
effects  of  a  single  load  block  or  processing  a 
preprogrammed  series  of  calculations; 

•  select  the  name  of  the  input  file  containing  the  load 
history; 

•  select  the  name  of  the  output  file; 

•  activate  and  deactivate  the  program's  sound. 

The  sound  function  mentioned  above  is  to  assist  the 
program's  user  when  processing  a  large  batch  of  calculations 
or  making  a  long  life  calculation.  When  activated,  it  will 
beep  every  time  the  program  finishes  processing  a  load  block 
and  sound  an  alarm  when  completely  finished.  From  the  second 
option  menu  the  user  has  the  choice  of  returning  to  the  first 
option  menu  by  pressing  ESCAPE  or  starting  the  program's 
calculation  subroutines  by  pressing  RETURN. 
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UP  .  Move  to  next  field 

DOWN  .  Move  to  previous  field 

LEPT/R1GMT  ....  Change  field  up/down 

F2  .  Turn  sound  to  ON 

ENTER  ....  Start  with  current  values 
ESCAPE  . .  Return  to  previous  screen 


Type  of  inputs  (stress  or  strain) 
Equation 

Fixed  /  Varing  n’  and  K' 
Calculation  Type 
Input  file  name 
Output  file  name 


(  Stress  ] 

[  Morrow’s  equation  ] 

[  Fixed  n‘  and  K’  ] 

[  Load  blocks  to  failure 
[  STRESS . DAT  ) 

[  OUTPUT . DAT  ] 


] 


Figure  2  Option  Menu  2 


The  other  menu  used  in  the  program  is  to  facilitate 
operator  entry  of  a  new  material.  Called  by  the  first  option 
selection  menu,  and  shown  in  Figure  3,  it  displays  and  allows 
the  operator  to  update  the  following  material  properties: 

•  ultimate  strength  (Su) 

•  yield  strength  (Sy) 

•  fatigue  yield  strength  (Sy' ) 

•  strength  coefficient  (K) 

•  cyclic  strength  coefficient  (K') 

•  strain  hardening  exponent  (n) 

•  cyclic  strain  hardening  exponent  (n') 

•  ductility  coefficient  (ef) 

•  fatigue  ductility  coefficient  (ef') 

•  strength  coefficient  (<7f) 
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•  fatigue  strength  coefficient  <crf' ) 

•  fatigue  strength  exponent  (b) 

•  fatigue  ductility  exponent  (c) 

•  endurance  strength  (Sf) 

•  modulus  of  elasticity  (E) 


UP  .  Hove  to  next 

field 

DOWN  .  Hove  to  previous 

field 

LEFT  or  RIGHT  ....  Enter  a  new 

value 

ENTER  ...  Return  with  current  values 

ESCAPE  ....  Quit  material  data 

entry 

Ultimate  Strength 

(Su  in  ksi ) 

[ 

0 

] 

Yield  Strength 

(Sy  in  ksi ) 

[ 

0 

J 

Fatigue  Yield  Strength 

(Sy ’  in  ksi ) 

[ 

0 

] 

Strength  Coeficient 

(K  in  ksi ) 

[ 

0 

) 

Cyclic  Strength  Coefficient 

(K ’  in  ksi ) 

[ 

0 

] 

Strain  Hardening  Exponent 

(n) 

[ 

0.00 

] 

Cyclic  Strain  Hardening  Exponent 

(n* ) 

[ 

0  00 

] 

Ductility  Coefficient 

(epsilon  f) 

[ 

0.00 

] 

Fatigue  Ouctility  Coefficient 

(epsi Ion  f ’) 

[ 

0.00 

] 

Strength  Coefficient 

(sigma  f  in  ksi ) 

t 

0 

] 

Fatigue  Strength  Coefficient 

(sigma  f ’  in  ksi ) 

[ 

0 

] 

Fatigue  strength  Exponent 

(b) 

[ 

0.00 

] 

Fatigue  Ductility  Exponent 

(c) 

[ 

0,00 

] 

Endurance  Strength 

(Sf  in  ksi ) 

[ 

0 

] 

Hodulus  of  Elasticity 

(E  in  ksi ) 

[ 

0 

] 

Figure  3  Material  Entry  Menu 


Upon  exiting  this  menu,  the  user  is  prompted  for  the  newly 
entered  material's  name  and  gives  the  option  of  saving  the  new 
material  in  the  material  data  base. 
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D.  LOAD  HISTORY  INPUT 

The  load  history  to  be  processed/analyzed  is  read  into  the 
program  from  the  operator  designated  input  file.  This  is 
accomplished  by  the  subroutine  LOADER.  Additionally,  this 
subroutine  determines  the  number  of  loads  in  the  sequence  and 
whether  the  first  increment  is  increasing  or  decreasing.  The 
load  input  file  is  simply  a  consecutive  list  of  either 
stresses  or  strains,  one  per  line,  in  a  DOS  compatible  format. 

E.  CALCULATIONS 

The  data  processing  takes  place  in  the  subroutine 
CRUNCHER.  In  this  subroutine  the  far- field  stress  or  strain 
load  sequence  is  turned  into  a  sequence  of  far- field  stress  or 
strain  changes,  according  to  what  was  initially  read  in.  If 
the  input  file  consisted  of  strain  loads,  the  monotonic 
stress-strain  equation,  (Equation  (6)  )  and  the  hysteresis 
equation,  (Equation  (7),)  is  used  to  convert  strain  changes  to 
stress  changes.  Next  using  Neuber's  rule,  (Equation  (5),)  the 
far- field  stress  changes  are  used  to  find  the  local  stress 
changes.  The  local  stress  changes  and  either  the  mean  or 
maximum  stress  are  then  used  in  one  of  the  three  available 
strain-life  equations,  (Equations  (2)  through  (4),)  to 
determine  the  number  of  cycles  to  failure  at  the  given  load 
levels.  This  life  is  then  used  in  conjunction  with  Miner's 
rule,  (Equation  (l))  to  determine  how  much  damage  has 
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accumulated,  or  how  much  of  the  fatigue  life  was  consumed,  and 
added  to  a  running  counter.  This  process  is  repeated  until 
all  the  loads  are  processed,  or  until  failure,  which  occurs 
when  the  damage  sums  to  one. 


F.  COMPUTATIONAL  SUBROUTINES 

Equations (2)  through  (6)  can  not  be  solved  directly  for 
the  needed  variables  and  require  an  iterative  numerical 
solution.  Newton's  method  of  successive  approximations  was 
chosen  as  the  simplest  and  most  efficient  method  for  this 
application.  Newton's  method  is  an  iterative  method  which 
requires  the  function  in  question  be  evaluated  at  an  initial 
estimate.  If  the  estimate  produces  excessive  error,  the 
function's  derivative  is  evaluated  at  the  estimated  value  and 
a  new  estimate  is  obtained  by  subtracting  from  the  old 
estimate  the  error  at  the  old  estimate  divided  by  the 
derivative  at  the  old  estimate.  [Reference  4] 

Subroutines  EQUATI0NS1,  EQUATIONS 2 ,  EQUATI0NS3 , 
EQUATI0NS4,  and  EQUATI0NS5  utilize  Newton's  method  to  solve 
equations  (2)  through  (5)  and  (7)  .  In  EQUATI0NS2 ,  EQUATI0NS3 , 
and  EQUATI0NS4 ,  which  solve  the  strain- life  equations, 
Newton's  method  was  modified  slightly  to  ensure  that  variable 
NNf  ^  being  solved  for  was  never  approximated  as  being 
negative.  (The  variable  Nnfif  which  is  raised  to  a  negative 
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power, 


would  produce  a  complex  number  which  Che  program 


cannoc  accept . ) 

To  allow  for  cycle  to  cycle  variation  in  K'  and  n'  the 
functions  xxkf  and  xxnf  were  created  to  redefine  the  values  of 
K'  and  n'  based  on  the  number  of  elapsed  cycles  (calculated  by 
the  function  xximf count)  . 

Since  QuickBasic  has  no  built  in  function  for  base  10 
logarithms,  the  function  LOGIO  was  incorporated  to  generate 
base  ten  logarithms  from  natural  logarithms  with  the 
relationship  shown: 


Logia(x)  = 


_  ln(x) 


In (10) 


(8) 


[Reference  5] . 


G.  DATA  OUTPUT 

Depending  on  the  processing  option  chosen,  the  output  will 
consist  of  either  a  step  by  step  print  out  of  each  variable's 
value  for  one  trip  through  the  input  load  block,  (Figure  4,) 
and  a  summary  of  the  chosen  options,  the  number  of  cycles  and 
load  blocks  until  failure,  or  just  the  output  summary,  (Figure 
5)  .  The  output  is  appended  to  the  output  file  so  that  any 
previous  data  in  the  file  will  not  be  lost. 

H.  CODE  VERIFICATION 


To  verify  the  programming  code,  the  FLP  was  run  for 
several  short  load  histories  (approximately  10  cycles) .  The 
first  test  sequences  were  constant  amplitude  loading  till 
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"Morrow's  equation" 

"  Fixed  n'  and  K'  " 

"Load  blocks  to  failure" 

"input  f ile: " , "testl" 

"output  file: ", "OUTPUT.DAT" 

"blocks :  " ,  7 

"i  counter 2482 

"reversal  count 35655 

"life  factor: ",1.000275410409951 


Figure  4  Sample  output  summary. 


index 

Stress 

del  t-aStrefls 

deltasiq 

siq 

siqO 

d“ltaepe 

Ntlf 

i 

45 

45 

45.00 

45.00 

22.50 

0.004369 

0 . 15Di09 

2 

20 

25 

25.00 

20.00 

32.50 

0  002428 

0  160*09 

3 

90 

70 

70.00 

90.00 

55.00 

0.006800 

0 . 66D+05 

4 

30 

60 

60.00 

30.00 

60.00 

0.005826 

0  14D*06 

5 

50 

20 

20.01 

50.00 

40.00 

0.001942 

0.180*09 

6 

10 

40 

40.01 

9.99 

30.00 

0.003884 

0.17D*09 

7 

40 

30 

30.02 

40.02 

25.00 

0.002915 

0 . 14D*09 

8 

90 

130 

130.00 

-89.98 

24.98 

0.012858 

0.18D*05 

9 

60 

150 

150.00 

60.02 

-14.98 

0.015193 

0.59D*04 

10 

-10 

70 

70.00 

-9.98 

25.02 

0 . 006800 

0 . 20D+06 

1 1 

90 

100 

100.00 

90.02 

40.02 

0.009748 

0 . 15D+05 

12 

40 

130 

130.00 

-39.98 

25 . 02 

0.012858 

0.58D*04 

13 

30 

70 

70.00 

30.02 

4.98 

0.006800 

0 . 150*09 

14 

0 

30 

30.02 

-0.01 

15.00 

0.002915 

0.15D*09 

15 

60 

60 

60 . 00 

59.99 

29.99 

0.005826 

0 . 12D+09 

IS 

10 

50 

50.01 

9.98 

34  99 

0.004856 

0 . 10D409 

17 

30 

20 

20 . 01 

29.99 

19.99 

0.001942 

0.110409 

18 

-30 

60 

60.00 

-30.01 

-0.01 

0.005826 

0.11D409 

19 

87 

117 

117.00 

86.99 

28.49 

0.011474 

0  91D*04 

20 

43 

44 

44 .00 

42.99 

64.99 

0.004272 

0.14D409 

21 

56 

13 

13.00 

55.99 

49.49 

0.001262 

0 . 1 1D+09 

22 

-76 

132 

132 . 00 

-76.01 

-10.01 

0.013078 

0 . 1 1D*05 

23 

4 

80 

80.00 

3.99 

-36.01 

0.007775 

0.15D409 

24 

5 

9 

9.01 

-5.02 

-0.52 

0.000875 

0.17D409 

25 

0 

5 

5.02 

-0.01 

-2.51 

0.000487 

0.20D409 

Figure  5  Sample  step  by  step  output. 


failure  and  later  load  sequences  were  variable  amplitude.  The 
identical  load  sequences  processed  by  hand  with  an  HP- 48  hand 
held  calculator  capable  of  solving  complex  equations 
numerically.  The  results  were  identical  to  four,  usually  f?vre 
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or  six,  significant  figures.  This  however  didn't  test  the 
programs  calculation  of  varying  strength  coefficient  and 
varying  strain  hardening  exponents.  To  verify  that  these  were 
being  correctly  calculated  the  program  was  run  in  the 
QuickBasic  environment  before  being  compiled  into  an 
executable,  (.EXE)  file.  When  the  program  is  being  run  in  the 
QuickBasic  environment,  the  operator  can  pause  execution  at 
any  time  and,  using  what's  referred  to  as  the  Immediate 
window,  execute  other  instructions  [Reference  6] .  Using  this 
feature  and  a  hand  held  calculator,  it  was  simple  to  pause  the 
program  at  random  intervals  between  1000  and  500,000  elapsed 
cycles  to  ensure  n'  and  K'  had  correct  values. 


III.  MATERIAL  DATA  BASE 


A.  CYCLIC  PROPERTIES 

In  most  metals  the  stress  -  strain  response  is  altered  by- 
repeated  loading.  Depending  on  the  material  and  it's  initial 
characteristics  (i.e.,  quenched  and  tempered,  or  annealed,)  a 
metal  might  either  cyclically  harden,  soften,  or  have  a  mixed 
behavior.  Typically  a  material  which  has  a  low  yield  strength 
relative  to  it's  ultimate  tensile  strength  will  undergo  strain 
hardening  while  a  hard  material  cyclically  softens.  This  is 
reflected  in  the  difference  between  a  material's  monotonic 
strength  coefficient  (K)  and  its  cyclic  strength  coefficient 
(K' )  as  well  as  the  strain  hardening  exponent  (n)  and  the 
cyclic  strain  hardening  exponent  (n' ) .  Figure  (6)  shows  the 
effects  of  cyclic  loading  on  copper. 

Normally,  when  performing  strain- life  calculations  K'  and 
n'  are  used  in  all  the  stress -strain  calculations.  For  low  to 
medium  cycle  fatigue,  (103  -  105  cycles,)  the  effect  of  the 
change  in  these  material  properties  is  undocumented.  To 
account  for  the  dynamic  nature  of  these  properties,  it  was 
postulated  that  they  could  be  expressed  as  a  function  of  the 
elapsed  cycles,  allowing  their  value  to  be  continually  updated 
from  cycle  to  cycle. 
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Figure  6  Hysteresis  response  of  copper.  (From  Ref.  1) 

B .  ARCHIVAL  DATA 

The  program  utilizes  a  separate  file  MAT. DAT  to  hold  the 
materials  available  to  the  program  without  manual  entry  by  the 
operator.  Appendix  z  shows  the  contents  of  the  file  MAT.  DAT. 

The  first  entry  in  the  material  data  base  indicates  the 
number  of  materials  in  the  original  data  base.  This  number 
couldn' t  be  changed  without  restructuring  the  data  base  into 
a  series  of  records  so  a  second  file  NEWCOUNT  was  created  just 
to  record  the  number  of  new  materials  that  were  added  to  the 
data  base.  The  number  of  materials  within  the  data  base  is 
needed  by  the  program  to  allow  the  option  menu  to  cycle 
through  all  the  materials. 

The  material  properties  used  in  the  data  base  came  from 
Metal  Fatigue  in  Engineering  by  H.  0.  Fuchs  [Reference  3]  . 


16 


The  data  base  contains  both  the  SI  and  British/American  values 
for  all  properties  and  when  new  materials  are  added  to  the 
data  base  both  the  SI  and  British/American  values  are 
recorded.  For  numerous  materials  Reference  3  doesn't  include 
values  for  K,  K'  and  n.  To  get  around  this  potential  problem, 
when  any  of  these  properties  are  missing  the  subroutine 
Loadmaterial ,  which  selects  the  proper  material  from  the  data 
base  and  initializes  the  variables  in  the  program  to  the 
appropriate  material  properties,  will  estimate  the  value  of 
the  missing  property  using  the  following  relations: 


K  = 


(9) 


[Reference  1]  . 

Review  of  the  selected  material's  properties  is 
accomplished  by  pressing  FI  while  the  first  option  menu  is 
displayed.  This  activates  the  display  used  for  entering  new 
materials,  which  will  display  the  properties  of  the  material 
that  is  already  loaded.  The  displayed  properties  can't  be 
modified  to  ensure  the  integrity  of  the  original  data  base. 
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IV.  LOAD  GENERATING  PROGRAM  LOAD GEN 


A.  LOAD  SPECTRUM  CONCEPT 

In  order  to  evaluate  the  effect  of  varying  n'  and  K'  on 
strain- life,  a  realistic  load  sequence  was  desired.  Readily 
available  was  the  load  spectrum  data  representing  the  g 
exceedences  of  a  0.9973  percentile  (three  sigma)  A- 6  aircraft, 
Table  I,  [Reference  2]. 


TABLE  1:  THREE  a  LOAD  SPECTRUM 


Data  of 

Data 

Corraclad 
FLt.  Hours 

40  PER 
1000.0 

SO 

1000.0 

SO 

1000.0 

7G 

1000.0 

12/31/89 

72  7 

165 

72 

0 

0 

1/31/90 

48  7 

186 

27 

0 

0 

2/28/90 

34  4 

74 

0 

0 

0 

3/31/90 

16  3 

74 

18 

0 

0 

4/30/90 

40  5 

234 

54 

0 

0 

5/31/90 

50.3 

90 

18 

0 

10 

6/30/90 

77  9 

32 

18 

0 

0 

7/31/90 

33  9 

42 

0 

0 

0 

8/31/90 

80  4 

127 

0 

0 

0 

9/30/90 

44  5 

32 

0 

0 

0 

10/31/90 

103  6 

90 

18 

0 

0 

11/30/90 

42  8 

111 

0 

0 

0 

12/31/90 

29  5 

27 

0 

0 

0 

7/31/91 

4  1 

106 

0 

0 

0 

8/31/91 

77  6 

58 

0 

0 

0 

9/30/91 

74  6 

42 

18 

0 

0 

10/31/91 

85  3 

260 

18 

0 

0 

1 1/30/91 

83.0 

228 

72 

47 

0 

Totals 

1000 

1978  0 

333  0 

47.5 

10.2 

The  number  of  exceedences  for  each  g  ~vei  was  entered 
into  the  program  LOADGEN.  This  program  .  -elect  a  random 
number  and  then  pick  a  corresponding  4g,  5g,  6g,  or  7g  load. 
The  random  number  generator  returned  a  number  between  0  and  1, 
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which  was  compared  to  the  percentage  of  g  occurrences  greater 
than  7,  then  greater  than  6,  and  so  on,  until  the  random 
number  is  less  than  or  equal  to  the  percentage  of  occurrences 
in  excess  of  a  certain  g  level.  The  algorithm  checks  to 
ensure  that  all  the  g  loads  at  that  level  haven't  been 
depleted  and  decrements  the  counter  for  the  respective  bin. 
If  all  the  loads  at  a  particular  g  level  are  depleted  and  a 
random  number  was  generated  picking  a  depleted  g  level,  then 
another  random  number  is  called  until  all  the  exceedences  have 
been  picked  from  the  four  g  levels. 

As  each  random  number  was  generated,  and  a  g  level 
exceedence  chosen,  the  program  immediately  calls  another 
random  number  to  determine  the  tenths  value  of  the  load  to  be 
inserted  into  the  sequence.  For  example,  if  a  5g  load  was 
picked  to  be  inserted  into  the  sequence,  the  second  random 
number  would  determine  if  it  were  a  value  of  5.0,  5.1,  5.2, 
etc.  up  to  5.9.  The  generated  values  were  then  sequenced  in 
an  output  file  for  temporary  storage.  With  the  assumption 
that  the  load  returned  to  1  g  at  the  end  of  every  cycle,  a 
sequence  like  the  one  in  Figure  7  was  generated. 

B.  A- 6  STRESS  SEQUENCE  GENERATION 

After  the  g  sequence  was  developed,  it  was  converted  to  a 
stress  sequence  for  processing  by  FLP.  Since  the  exact 
correlation  between  g  load  and  the  stress  it  produces  on  the 
A- 6  wasn't  known,  the  standard  Navy  design  criteria  was  used 
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to  calculate  the  highest  acceptable  corresponding  stress 
values.  The  Navy's  design  criteria  is  such  that  the  yield 
stress  can't  be  exceeded  at  the  design  g  limit,  and  the 
ultimate  stress  can't  be  exceeded  at  1.5  times  the  design  g 
limit.  Since  the  design  g  limit  of  the  A-6  is  +6.5,  in  the 
worst  case  either  a  g  loading  of  6.5  corresponds  to  the 
material's  yield  strength,  or  a  g  loading  of  9.75  corresponds 
to  the  material's  ultimate  strength.  Given  both  the  yield  and 
ultimate  strengths  of  the  material,  the  program  selects  the 
limiting  case  and  uses  the  corresponding  g  to  stress  ratio  to 
convert  the  series  of  g  loadings  into  a  series  of  stress 
loadings.  Since  the  g  count  information  only  gives  maximums 
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and  no  minimums,  the  load  profile  was  assumed  to  return  to  1 
g  between  each  of  the  peaks. 

As  the  file  of  g  loads  was  transformed  into  stress  loads, 
the  stresses  were  written  sequentially  into  a  file  for  use  by 
FLP. 
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V. 


APPROXIMATION  OF  K'  AND  n 


A.  EXPERIMENTAL  EFFORTS 

To  account  for  the  cyclic  change  in  K'  and  n'  ,  these 
values  were  experimentally  determined  in  test  specimens  whose 
life  under  a  given  cyclic  strain  load  had  already  been 
determined.  Specimens  made  from  aluminum  7075 -T6,  shown  in 
Figure  8,  were  tested  in  cyclic  strain  by  LT  Byron  Smith  in  an 
investigation  of  the  effects  of  means  strain  on  strain  life 
[reference  7].  Under  fully  reversed  strain  loading  of  ±.007 
in/in,  test  specimens  were  cycled  to  10%,  20%,  30%,  and  40%  of 
their  predetermined  total  life.  These  specimens,  and  two 
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specimens  without  any  cyclic  preloading,  were  then  subjected 
to  a  uniaxial  stress -strain  test  to  determine  the  properties 
in  question. 

The  testing  was  performed  on  the  Naval  Postgraduate 
School,  Mechanical  Engineering  Department's  MTS  machine.  Data 
correlating  to  load  and  displacement  were  recorded 
continuously  throughout  the  tests.  This  data  was  reduced  to 
stress  and  strain,  plotted  in  Figure  9,  and  then  processed  to 
determine  the  stress  coefficients  and  strain  hardening 
exponents  at  the  time  of  testing. 

Two  methods  of  determining  the  properties  in  question  were 
attempted.  In  the  first  the  strain  hardening  exponents  (n) 
were  calculated  using  the  following  relationship  from 
reference (1) : 

at  necking:  n  =  ln(l  +  e)  (10) 

where  e  is  the  engineering  strain. 

Knowing  n,  K  can  be  determined  with  the  relationship: 

« -*■<«»>*  -  an 

where  a  is  the  true  stress  and  ep  is  the  true  elastic  strain. 
Taking  into  consideration  the  following  relationships: 

e  =e  -e  (12) 

cp  ^ total 


£ 


E 


(13) 
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(14) 


e£otai  =  ln(i  +  e) 

a  =  5(1  +  e) ;  (15) 

where  eg  is  the  true  elastic  strain,  and  S  is  the  engineering 
stress,  this  expression  can  be  formed: 

(16) 

The  resulting  values  for  n'  and  K'  appear  in  Table  2. 


|ln(l  +  e)  - 


TABLE  2:  DATA  REDUCTION  RESULTS  -  METHOD  1 


n 

K  (ksi) 

0%  prestrain  #1 

.160 

131 

0%  prestrain  #2 

.155 

130 

10%  prestrain 

.020 

93 

20%  prestrain 

.155 

131 

30%  prestrain 

.  151 

130 

40%  prestrain 

.155 

132 
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The  experimentally  derived  stress  coefficient  and  strain 
hardening  data  doesn't  lend  itself  to  generating  a  function 
based  on  elapsed  cycles  or  used  up  strain- life.  Looking  at 
more  of  LT  Smith's  work  with  strain  fatigue  in  aluminum  7075- 
T6 ,  shows  a  large  deviation  in  the  test  results  for  strain 
life,  (Figure  10)  .  Because  of  the  large  deviation  in  the  life 
of  a  specimen,  there  are  large  deviations  in  specified 
percentages  of  the  sample's  life.  This  indicates  the  need  to 
test  a  large  number  of  samples  to  get  results  for  K  and  n  of 
any  statistical  significance. 


Figure  9  Stress -Strain  Curves 
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LOGNORMAL  PROBABILfTY  PLOT,  N=20 


1000  2000  4000  6000  POOO 


Figure  10  Strain- Life  Deviation 


B.  ESTIMATION  OF  K'  AND  n' 

Lacking  statistically  significant  test  results,  for 
purposes  of  exercising  this  program  both  the  strength 
coefficient  and  the  strain  hardening  exponent  were  postulated 
to  be: 

•  the  monotonic  values  when  less  than  1000  cycles  had 
elapsed, 

•  the  cyclic  values  when  more  than  500,000  cycles  had 
elapsed,  (half  the  endurance  life,)  and 

•  a  log-normal  function  of  the  number  of  elapsed  cycles  when 
between  1,000  and  500,000. 


To  calculate  these  properties  between  1000  and  500,000  cycles 
the  following  equations  were  used: 

)L0?-° (JC>  -  “W1000*)  - « 

(17) 

Kvazi  =  (  Log10  (5000ofy  -  Log10  (1000)  \LOg™  {x)  ~  ^SioUOOO))  +  K 

(18) 

For  aluminum  7075 -T6  these  equations  are  plotted  in  Figure  11. 
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Figure  11  K'  and  n’  Functions 
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VI .  COMPUTATIONS 


A.  COMPARISON  PROCEDURE 

To  compare  the  strain- lives  calculated  by  using  the  cyclic 
properties  with  those  obtained  by  varying  the  properties,  a 
set  of  input  load  files  were  generated  using  the  three  a 
limits  of  g  load  data  for  the  A- 6  and  the  LOADGEN  program. 
Each  file  contained  approximately  4800  points,  representing 
1000  hours  of  flight  time.  The  subroutine  BATCH  was  setup  to 
repeat  the  strain- life  calculation  for  each  of  the  four  load 
files,  using  each  of  the  three  strain- life  equations,  and  the 
four  following  calculation  methods: 

•  fixed  n'  and  fixed  K' 

•  fixed  n'  and  variable  K' 

•  variable  n'  and  fixed  K' 

•  variable  n'  and  variable  K' 

This  amounted  to  48  sets  of  computations,  the  results  of  which 
are  shown  in  Table  3 . 

B.  DISCUSSION  OF  RESULTS 

The  data  in  Table  3  show  a  reduction  in  the  calculated 
fatigue  life  of  15  to  30%  when  both  the  strength  coefficient 
and  the  strain  hardening  exponent  were  varied  from  their 
monotonic  values  to  their  cyclic  values.  It  is  also  clear 
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-hat  the  strength  coefficient  has  a  much  greater  impact  on  the 
fatigue  life  then  does  the  strain  hardening  exponent. 
Variation  of  the  strength  coefficient  consistently  produced  a 
significant  drop  in  the  fatigue  life,  while  varying  the  strain 
hardening  exponent  very  slightly  raised  the  calculated  fatigue 
life . 

All  three  of  the  strain- life  equations  used  to  make  the 
calculations  behaved  the  same,  though  it  seemed  that  the 
respective  importance  of  the  strain  harding  exponent  relative 
to  the  strength  coefficient  varied  between  the  equations. 

The  effect  of  varying  the  strength  coefficient  can  be 
anticipated.  From  Equation(6),  a  reduced  value  of  K' 
translates  into  greater  strain  values  locally,  and  greater 
strain  means  reduced  life.  However  reduced  n'  values  also 
mean  larger  strains  and  consequent  reductions  in  life. 

The  effect  of  varying  only  the  strain  hardening  could  also 
be  anticipated,  but  it's  impact,  compared  to  that  of  the 
strength  coefficient  was  unexpected. 

It  should  be  noted  that  these  calculations  assumed 
constant  values  for  b  and  c  in  Equations  (2)  through  (4)  , 
where  in  reality  they  would  change.  However  this  example  has 
allowed  this  unique  segment  of  the  program  to  be  exercised  and 
to  gain  some  preliminary  estimates  of  the  effect  varying  these 
parameters  has  on  fatigue  life. 
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TABLE  3:  CYCLES  TO  FAILURE 


Equation 

P  rn rpssi n  g 

Input 

file : 

Option 

TESTO 

TEST1 

TEST2 

TEST3 

Morrow' s 

Fix  n' ;  Fix  K' 

50554 

49570 

49708 

50221 

Equation 

Var  n' ;  Var  K' 

38723 

33665 

33661 

37899 

.... 

Var  n' ;  Fix  K' 

50783 

" 

49737 

50000 

50453 

Fix  n' ;  Var  K' 

38723 

33665 

33661 

38169 

Smith- 

Fix  n' ;  Fix  K' 

39018 

38404 

38543 

37899 

Watson- 

Topper 

Var  n' ;  Var  K' 

33985 

29979 

31340 

33114 

Var  n' ;  Fix  K' 

39162 

38635 

38779 

38847 

Fix  n' ;  Var  K' 

33933 

28927 

30310 

32072 

Mans on- 

Fix  n' ;  Fix  K' 

47428 

46483 

46838 

47346 

Halford 

Var  n' ;  Var  K' 

34516 

32166 

33661 

33729 

Var  n' ;  Fix  K' 

47587 

46924 

47320 

47429 

Fix  n' ;  Var  K' 

33985 

28927 

33661 

33729 
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VII.  DISCUSSION/CONCLUSIONS 


The  FLP  originally  conceived  has  been  wiitten  and  appears 
to  function  well.  It  possesses  flexibity  in  the  manner  in 
which  it  performs  fatigue  life  calculations,  using  any  one  of 
three  documented  strain- life  equations  which  account  for  mean 
stress  effects.  The  initial  exploration  done  with  it  supports 
tne  concept  that  fatigue  life,  especially  in  the  lower  end  of 
the  fatigue  spectrum,  might  be  better  calculated  if  changes  in 
material  properties  are  accounted  for.  To  do  this,  however, 
material  properties  need  to  be  well  defined  as  functions  of 
elapsed  cycles.  The  algorithms  used  in  the  FLP  varied  only 
two  of  the  material's  properties  as  a  function  of  elapsed 
cycles,  K'  and  n'  .  Other  material  properties  could  be  varied 
in  addition  to  the  ones  currently  varied  by  the  FPL.  Four 
candidates  are:  the  fracture  strength  coefficient  (af) ,  the 
ductility  coefficient  (ef)  ,  the  fatigue  strength  exponent  (b) 
and  the  fatigue  ductility  exponent  (c)  .  These  coefficients 
have  been  shown  to  vary  with  elapsed  cycles  as  well,  and  their 
effect  should  be  explored  in  future  tests  and  incorporated 
into  FLP. 
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APPENDIX  A  FLP  PROGRAM  LISTING 


A.  MAIN  PROGRAM 

DEFSNG  A-Z 
DEFINT  I-P 
DECLARE  SUB  BATCH  () 

DECLARE  SUB  LOADER  () 

DECLARE  SUB  CRUNCHER  (} 

DECLARE  SUB  EQUATIONS1  () 

DECLARE  SUB  EQUATIONS 2  () 

DECLARE  SUB  EQUATIONS3  0 
DECLARE  SUB  EQUATIONS4  0 
DECLARE  SUB  EQUATIONS 5  0 
DECLARE  SUB  OPTMENU1  () 

DECLARE  SUB  OPTMENU2  () 

DECLARE  SUB  MATMENU  () 

DECLARE  SUB  UPDATEMENU  0 
DECLARE  SUB  UPDATEMENU2  {) 

DECLARE  SUB  Loadmaterial  (index) 

DECLARE  SUB  NEWMAT  (Ky) 

DECLARE  SUB  OUTPUTER  () 

DECLARE  SUB  HEADER  () 

DECLARE  SUB  DATADUMP  () 

DECLARE  SUB  PetersonKf  () 

DECLARE  SUB  NeuberKf  () 

DECLARE  SUB  GetConfig  () 

DECLARE  SUB  TYPIN  (cstring$,  valu) 

DECLARE  SUB  TYPINSTRING  (cstring$,  stringname$ ) 

DECLARE  SUB  Klaxon  (Hi%,  low%) 

DECLARE  FUNCTION  xxnf !  () 

DECLARE  FUNCTION  xxKf !  () 

DECLARE  FUNCTION  xxNNf count &  () 

DECLARE  FUNCTION  Rotated  (Lower,  Upper  AS  INTEGER,  present,  Inc) 
DECLARE  FUNCTION  LOGIO !  (value!) 


Main  Program 


REM  $ DYNAMIC 

'  Constants  for  best  available  screen  mode 

CONST  VGA  =12 

CONST  MCGA  =13 

CONST  EGA256  =  9 

CONST  EGA64  =  8 

CONST  MONO  =10 

CONST  HERC  =  3 

CONST  CGA  =  1 


'  User-defined  type  to  hold  information  about  the  mode 
TYPE  Config 

Scrn  AS  INTEGER 

Colors  AS  INTEGER 

Atribs  AS  INTEGER 

XPix  AS  INTEGER 

YPix  AS  INTEGER 

TOOL  AS  INTEGER 

TROW  AS  INTEGER 

END  TYPE 

DIM  VC  AS  Config 


'  variables  shared  throughout  the  program 

COMMON  SHARED  YD,  i,  top,  Kt  AS  SINGLE,  Ktf  AS  SINGLE,  option4  AS  STRING 
COMMON  SHARED  matcount,  matindex,  matname  AS  STRING,  Su,  Sy,  Syf,  K  AS  SINGLE 
COMMON  SHARED  Kf  AS  SINGLE,  n  AS  SINGLE,  nf  AS  SINGLE,  equationname  AS  STRING 
COMMON  SHARED  epsf,  epsff,  sigf,  sigff,  b,  c,  Sf,  SfSu,  E,  newstuff 
COMMON  SHARED  flagl  AS  INTEGER,  flag2  AS  INTEGER,  flag3  AS  INTEGER 
COMMON  SHARED  flag4  AS  INTEGER,  flag5  AS  STRING,  batchno  AS  INTEGER 
COMMON  SHARED  nKtype  AS  STRING,  calctype  AS  STRING,  IgNNf count  AS  SINGLE 
COMMON  SHARED  usedlife  AS  DOUBLE,  block  AS  INTEGER,  lasti  AS  INTEGER 
COMMON  SHARED  NNf count  AS  LONG,  xKf  AS  SINGLE,  xnf  AS  SINGLE 
COMMON  SHARED  inputfile  AS  STRING,  outputfile  AS  STRING,  sigmax  AS  SINGLE 
COMMON  SHARED  laststress  AS  SINGLE,  strainmax  AS  SINGLE,  lasteps  AS  SINGLE 
COMMON  SHARED  lastsig  AS  SINGLE 

size  =  5000 

DIM  SHARED  stress(size)  AS  SINGLE,  strain (size) ,  sig(size)  AS  SINGLE 
DIM  SHARED  deltaeps (size)  AS  SINGLE,  NNf (size)  AS  DOUBLE,  sigO(size) 

DIM  SHARED  deltaStress (size)  AS  SINGLE,  deltastrain (size)  AS  SINGLE 
DIM  SHARED  deltasig ( size ) ,  eps(size) 


'  User-defined  type  to  hold  information  about  the  selected  options  (menul) 

TYPE  Optionsl 

units  AS  STRING  *  4 
material  AS  INTEGER 
Kt  AS  SINGLE 

Ktf  AS  INTEGER 

option5  AS  INTEGER 

END  TYPE 

DIM  Menu  AS  Optionsl 

'  User-defined  type  to  hold  information  about  the  selected  options  (menu2) 

TYPE  Options2 

inputs  AS  STRING  *  6 
equation  AS  INTEGER 

option3  AS  INTEGER  'nKtype  selection 

option4  AS  INTEGER  'calculation  type 

END  TYPE 

DIM  menu2  AS  Options2 

'  Error  variables  to  check  screen  type 

DIM  InitRows  AS  INTEGER,  BestMode  AS  INTEGER,  Available  AS  STRING 
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'  Initialize  variables 

Menu. units  =  "Brit" 

Menu .material  =  21 

Menu . Kt  =  l 

Menu.Ktf  =  l 

Menu. options  =  1 

option4  =  "Manually  Entered  " 

Kt  =  1 
Ktf  =  0 

menu2 . inputs  =  "Stress" 
menu2 . equation  =  l 
menu2 ,option3  =  1 
menu2 ,option4  =  1 

equationname  =  "Morrow' s  equation" 
nKtype  =  "  Fixed  n'  and  K'  " 
calctype  =  "Load  blocks  to  failure" 
inputfile  =  "STRESS.DAT" 
outputfile  =  " OUTPUT . DAT " 


flag2  =  l 
flag3  =  0 
flag4  =  0 

flags  =  "OFF"  '  initialises  sound  on 

batchno  =  0 

sigmax  =  0 

strainmax  =  0 

laststress  =  0 

lastsig  =  0 

lasteps  =  0 

'  Get  best  configuration  and  set  initial  graphics  mode  to  it 

GetConf ig 

VC.Scrn  =  BestMode 


DO  '  endless  do  loop  for  main  program  (exited  from  menu  l) 

'  Selecti  n  of  material/constants 

IF  flag4  =  0  THEN 
DO 

CALL  0PTMENU1 
CALL  0PTMENU2 
LOOP  UNTIL  f lag2  =  1 
END  IF 

IF  menu2 ,option4  =  3  THEN 
batchno  =  batchno  +  1 
BATCH 
END  IF 

CALL  LOADER 
CALL  CRUNCHER 

CALL  OUTPUTER 
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IF  flag4  =  0  THEN 

IF  flag5  =  "OFF"  THEN  PRINT  "Finished  Processing . press  any  key" 

IF  flags  =  "OFF"  THEN  CALL  Klaxon (987,  329) 

CLS 

ELSE 

PRINT  USING  "Finished  Processing  Batch  ##";  batchno 
END  IF 

sigmax  =  0 
strainmax  =  0 
laststress  =  0 
lasteps  =  0 
lastsig  =  0 

LOOP  UNTIL  cowscomehome 
END 


'  Error  trap  to  make  program  screen  independent 
VideoErr : 

SELECT  CASE  BestMode  '  Fall  through  until  something  works 
CASE  VGA 

BestMode  =  MCGA 
Available  =  "12BD" 

CASE  MCGA 

BestMode  =  EGA25G 
Available  =  "12789" 

CASE  EGA256 

BestMode  =  CGA 
Available  =  "12" 

CASE  CGA 

BestMode  =  MONO 
Available  =  "A" 

CASE  MONO 

BestMode  =  HERC 
Available  =  "3" 

CASE  ELSE 

PRINT  "Sorry.  Graphics  not  available.  " 

END 

END  SELECT 

RESUME 


'  Trap  to  detect  64K  EGA 
EGAErr : 

BestMode  =  EGA64 
Available  =  "12789" 
RESUME  NEXT 


'  Trap  to  detect  insufficient  memory 
MemErr : 

LOCATE  22,  1 

PRINT  "Out  of  memory" 

RESUME  NEXT 
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'  Trap  to  determine  initial  number  of  rows  so  they  can  be  restored 
RowErr : 

IF  Ini t Rows  =50  THEN 
InitRows  =  43 
RESUME 
ELSE 

InitRows  =  25 
RESUME  NEXT 
END  IF 


=  end  of  main  program  +  error  trapping  ======== 


B .  SUBROUTINE  BATCH 


REM  $STATIC 

'  Subroutine  to  access  a  file  for  repeated  runs 
'  (currently  set  to  process  4  file,  with  3  eqns,  4  methods 


SUB  BATCH 

SHARED  menu2  AS  Options2 

IF  batchno  =  1  THEN 

inputfile  =  "testaa" 
nKtype  =  "  Fixed  n'  and  K'  " 
menu2 .option3  =  1 
menu2 .equation  =  1 

equationname  =  "Morrow' s  Equation" 
ELSE IF  batchno  =  2  THEN 
inputfile  =  "testaa" 
nKtype  =  "Variable  n'  and  K' " 
menu2 . option3  =  2 
menu2 .equation  =  1 

equationname  =  "Morrow' s  Equation" 
ELSE IF  batchno  =  3  THEN 
inputfile  =  "testaa" 
nKtype  =  "Variable  n'  and  fixed  K'  " 
menu2 . option3  =  3 
menu2 .equation  =  1 
equationname  =  "Morrow's  Equation" 
ELSEIF  batchno  =  4  THEN 
inputfile  =  "testaa" 
nKtype  =  "Fixed  n'  and  variable  K' " 
menu2 . option3  =  4 
menu2 .equation  =  1 
equationname  =  "Morrow's  Equation" 
ELSEIF  batchno  =  5  THEN 
inputfile  =  "testbb" 
nKtype  =  "  Fixed  n'  and  K'  " 
menu2 .option3  =  1 
menu2 .equation  =  1 
equationname  =  "Morrow's  Equation" 
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ELSEIF  batchno  =  6  THEN 
inputfile  =  "testbb" 
nKtype  =  "Variable  n'  and  K'  " 
menu2 . option3  =  2 
menu2 . equation  =  1 
equationname  =  "Morrow's  Equation" 
ELSEIF  batchno  =  7  THEN 
inputfile  =  "testbb" 
nKtype  =  "Variable  n'  and  fixed  K'  " 
menu2 . option3  =  3 
menu2 . equation  =  1 
equationname  =  "Morrow's  Equation" 
ELSEIF  batchno  =  8  THEN 
inputfile  =  "testbb" 
nKtype  =  "Fixed  n'  and  variable  K' " 
menu2 ,option3  =  4 
menu2 . equation  =  1 
equationname  =  "Morrow's  Equation" 
ELSEIF  batchno  =  9  THEN 
inputfile  =  "testcc" 
nKtype  =  "  Fixed  n'  and  K'  " 
menu2 ,option3  =  1 
menu2 .equation  =  l 
equationname  =  "Morrow's  Equation" 
ELSEIF  batchno  =  10  THEN 
inputfile  =  "testcc" 
nKtype  =  "Variable  n'  and  K' " 
menu2 ,option3  =  2 
menu2 .equation  =  1 
equationname  =  "Morrow' s  Equation" 
ELSEIF  batchno  =  11  THEN 
inputfile  =  "testcc" 
nKtype  =  "Variable  n'  and  fixed  K’  " 
menu2 ,option3  =  3 
menu2 . equation  =  1 
equationname  =  "Morrow’ s  Equation" 
ELSEIF  batchno  =  12  THEN 
inputfile  =  "testcc" 
nKtype  =  "Fixed  n'  and  variable  K' " 
menu2 ,option3  =  4 
menu2 . equation  =  l 
equationname  =  "Morrow' s  Equation" 
ELSEIF  batchno  =  13  THEN 
inputfile  =  "testdd" 
nKtype  =  "  Fixed  n'  and  K'  " 
menu2 ,option3  =  1 
menu2 .equation  =  l 
equationname  =  "Morrow' s  Equation" 
ELSEIF  batchno  =  14  THEN 
inputfile  =  "testdd" 
nKtype  =  "Variable  n'  and  K'  " 
menu2 ,option3  =  2 
menu2 .equation  =  l 
equationname  =  "Morrow' s  Equation" 
ELSEIF  batchno  =  15  THEN 
inputfile  =  "testdd" 
nKtype  =  "Variable  n'  and  fixed  K'  " 
menu2 ,option3  =  3 
menu2 .equation  =  l 
equationname  =  "Morrow' s  Equation" 
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ELSEIF  batchno  =  16  THEN 
inputfile  =  "testdd" 
nKtype  =  "Fixed  n'  and  variable  K' " 
menu2 .option3  =  4 
menu2 . equation  =  l 
equationname  =  "Morrow's  Equation" 
ELSEIF  batchno  =  17  THEN 
inputfile  =  "testaa" 
nKtype  =  "  Fixed  n'  and  K'  " 
menu2 . option3  =  1 
mcnu2 . equation  =  2 

equationname  =  "Smith-Watson-Topper" 
ELSEIF  batchno  =  18  THEN 
inputfile  =  "testaa” 
nKtype  =  "Variable  n'  and  K' " 
menu2 .option3  =  2 
menu2 .equation  =  2 

equationname  =  "Smith-Watson-Topper" 
ELSEIF  batchno  =  19  THEN 
inputfile  =  "testaa" 
nKtype  =  "Variable  n'  and  fixed  K'  " 
menu2 ,option3  =  3 
menu2 .equation  =  2 

equationname  =  "Smith-Watson-Topper" 
ELSEIF  batchno  =  20  THEN 
inputfile  =  "testaa" 
nKtype  =  "Fixed  n'  and  variable  K' " 
menu2 .option3  =  4 
menu2 . equation  =  2 

equationname  =  "Smith-Watson-Topper" 
ELSEIF  batchno  =21  THEN 
inputfile  =  "testbb" 
nKtype  =  "  Fixed  n'  and  K'  " 
menu2 . option3  =  1 
menu2 .equation  =  2 

equationname  =  "Smith-Watson-Topper" 
ELSEIF  batchno  =  22  THEN 
inputfile  =  "testbb" 
nKtype  =  "Variable  n'  and  K' " 
menu2 .option3  =  2 
menu2 .equation  =  2 

equationname  =  "Smith-Watson-Topper” 
ELSEIF  batchno  =23  THEN 
inputfile  =  "testbb" 
nKtype  =  "Variable  n'  and  fixed  K'  " 
menu2 ,option3  =  3 
menu2 .equation  =  2 

equationname  =  "Smith-Watson-Topper" 
ELSEIF  batchno  =  24  THEN 
inputfile  =  "testbb" 
nKtype  =  "Fixed  n'  and  variable  K' " 
menu2 .option3  =  4 
menu2 .equation  =  2 

equationname  =  "Smith-Watson-Topper" 
ELSEIF  batchno  =  25  THEN 
inputfile  =  "testcc" 
nKtype  =  "  Fixed  n'  and  K'  " 
menu2 .option3  =  1 
menu2 .equation  =  2 

equationname  =  "Smith-Watson-Topper" 
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ELSEIF  batchno  =  26  THEN 
inputfile  =  "testcc" 
nKtype  =  "Variable  n'  and  K' " 
menu2 . opcion3  =  2 
menu2 .equation  =  2 

equationname  =  "Smith-Watson-Topper" 
ELSEIF  batchno  =  27  THEN 
inputfile  =  "testcc" 

nKtype  =  "Variable  n'  and  fixed  K'  " 
menu2 .option3  =  3 
menu2 . equation  =  2 

equationname  =  "Smith-Watson-Topper" 
ELSEIF  batchno  =28  THEN 
inputfile  =  "testcc" 
nKtype  =  "Fixed  n'  and  variable  K' " 
menu2 . option3  =  4 
menu2 .equation  =  2 

equationname  =  "Smith-Watson-Topper" 
ELSEIF  batchno  =  29  THEN 
inputfile  =  "testdd" 
nKtype  =  "  Fixed  n'  and  K'  " 
menu2 . option3  =  2 
menu2 .equation  =  1 

equationname  =  "Smith-Watson-Topper" 
ELSEIF  batchno  =30  THEN 
inputfile  =  "testdd" 
nKtype  =  "Variable  n'  and  K' " 
menu2 .option3  =  2 
menu2 .equation  =  2 

equationname  =  "Smith-Watson-Topper" 
ELSEIF  batchno  =31  THEN 
inputfile  =  "testdd" 
nKtype  =  "Variable  n'  and  fixed  K'  " 
menu2 . option3  =  3 
menu2 .equation  =  2 

equationname  =  "Smith-Watson-Topper" 
ELSEIF  batchno  =  32  THEN 
inputfile  =  "testdd" 
nKtype  =  "Fixed  n'  and  variable  K' " 
menu2 . option3  =  4 
menu2 .equation  =  2 

equationname  =  "Smith-Watson-Topper" 
ELSEIF  batchno  =33  THEN 
inputfile  =  "testaa" 
nKtype  =  "  Fixed  n'  and  K'  " 
menu2 ,option3  =  1 
menu2 .equation  =  3 
equationname  =  "Manson-Halford" 
ELSEIF  batchno  =  34  THEN 
inputfile  =  "testaa" 
nKtype  =  "Variable  n'  and  K' " 
menu2 .option3  =  2 
menu2 .equation  =  3 
equationname  =  "Manson-Halford" 
ELSEIF  batchno  =35  THEN 
inputfile  =  "testaa" 
nKtype  =  "Variable  n'  and  fixed  K'  " 
menu2 .option3  =  3 
menu2 .equation  =  3 
equationname  =  "Manson-Halford" 
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ELSE IF  batchno  =36  THEN 
inputfile  =  "testaa" 
nKtype  =  "Fixed  n'  and  variable  K' " 
menu2 .option3  =  4 
menu2 .equation  =  3 
equationname  =  "Manson-Half ord" 
ELSEIF  batchno  =  37  THEN 
inputfile  =  "testbb" 
nKtype  =  "  Fixed  n'  and  K'  " 
menu2 ,option3  =  1 
menu2 .equation  =  3 
equationname  -  "Manson-Halford" 
ELSEIF  batchno  =38  THEN 
inputfile  =  "testbb" 
nKtype  =  "Variable  n'  and  K' " 
menu2 . option3  =  2 
menu2 .equation  =  3 
equationname  =  ” Manson-Half ora" 
ELSEIF  batchno  =39  THEN 
inputfile  =  "testbb" 
nKtype  =  "Variable  n'  and  fixed  K'  ” 
menu2 . option3  =  3 
menu2 . equation  =  3 
equationname  =  "Manson-Halford" 
ELSEIF  batchno  =40  THEN 
inputfile  =  "testbb" 
nKtype  =  "Fixed  n'  and  variable  K' " 
menu2 .option3  =  4 
menu2 .equation  =  3 
equationname  =  "Manson-Halford" 
ELSEIF  batchno  =  41  THEN 
inputfile  =  "testcc" 
nKtype  =  "  Fixed  n'  and  K'  " 
menu2 . option3  =  1 
menu2 .equation  =  3 
equationname  =  "Manson-Halford" 
ELSEIF  batchno  =42  THEN 
inputfile  =  "testcc" 
nKtype  =  "Variable  n'  and  K' " 
menu2 .option3  =  2 
menu2 .equation  =  3 
equationname  =  "Manson-Halford" 
ELSEIF  batchno  =  43  THEN 
inputfile  =  "testcc" 
nKtype  =  "Variable  n'  and  fixed  K'  " 
menu2 . option3  =  3 
menu2 . equation  =  3 
equationname  =  "Manson-Halford" 
ELSEIF  batchno  =  44  THEN 
inputfile  =  "testcc" 
nKtype  =  "Fixed  n'  and  variable  K' " 
menu2 . option3  =  4 
menu2 .equation  =  3 
equationname  =  "Manson-Halford" 
ELSEIF  batchno  =45  THEN 
inputfile  =  "testdd" 
nKtype  =  "  Fixed  n'  and  K'  " 
menu2 ,option3  =  l 
menu2 .equation  =  3 
equationname  =  "Manson-Halford" 
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ELSEIF  batchno  =  46  THEN 
inputfile  =  "testdd" 
nKtype  =  "Variable  n'  and  K' " 
menu2 . option3  =  2 
menu2 . equation  =  3 
equationname  =  "Manson-Halford" 
ELSEIF  batchno  =47  THEN 
inputfile  =  "testdd" 
nKtype  =  "Variable  n'  and  fixed  K'  " 
menu2 .option3  =  3 
menu2 . equation  =  3 
equationname  =  "Manson-Halford" 
ELSEIF  batchno  =48  THEN 
inputfile  =  "testdd" 
nKtype  =  "Fixed  n'  and  variable  K' " 
menu2 .option3  =  4 
menu2 . equation  =  3 
equationname  =  "Manson-Halford" 
flag4  =  0 
END  IF 
END  SUB 


C.  SUBROUTINE  CRUNCHER 


REM  $DYNAMIC 

'  subroutine  to  calculate  Strian-Life  using  the  various  methods 

SUB  CRUNCHER 

DIM  lastload  AS  SINGLE 
SHARED  menu2  AS  0ptions2 

OPEN  outputfile  FOR  APPEND  AS  #7 

'initialize  variables 
usedlife  =  0 
block  =  -1 
lastload  =  0 

DO  '  loop  until  life  is  used  up  if  single  block  option  not  chosen 

block  =  block  +  l 

IF  block  =  0  THEN  'calculate  delta  loads  and  print  output  header  if 

appropriate 

IF  menu2 . option4  =  2  THEN  HEADER  '  output  headers 

IF  menu2. inputs  =  "Stress"  THEN 

'  calculation  of  farfield  deltastress  array 

FOR  ii  =  1  TO  top 

deltaStress { ii )  =  ABS (stress (ii)  -  lastload) 
lastload  =  stress (ii) 

NEXT  ii 

ELSE 
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'  calculation  of  far  field  deltastrain  array 

FOR  ii  =  1  TO  top 

deltastrain (ii)  =  ABS (strain ( ii )  -  lastload) 

lastload  =  strain (ii) 

NEXT  ii 
END  IF 


END  IF 


FOR  i 


load  blocks  processed 


1  TO  top 

NNfcount  = 
IgNNf count 
xKf  =  xxKf 
xnf  =  xxnf 


'  set  properties  based  on  elapsed  cycles 
xxNNf count 

=  LOGIO (CSNG( NNfcount ) ) 


'  turn  farfield  strains  into  stresses 
IF  menu2. inputs  =  "Strain"  THEN  CALL  EQUATIONS 5 


find  local  Stress  associated  with  farfield  stress 


CALL  EQUATIONS 1 

'  calcalate  mean  local  stress 
sig(i)  =  lastsig  +  ((-1)  A  (flagl  +  i) )  *  deltasig(i) 
sigO(i)  =  (lastsig  +  sig(i))  /  2 
lastsig  =  sig(i) 


'  finding  epsilon  and  delta  epsilon 
IF  sig(i)  >  sigmax  THEN 

'  mopnotonic  equation  to  find  local  strains 
eps(i)  =  sig(i)  /  E  +  (sig(i)  /  (xKf ) )  (1  /  xnf) 

deltaeps(i)  =  ABS(eps(i)  -  lasteps) 
lasteps  =  eps(i) 
sigmax  =  sig(i) 

ELSE 


'  hysteresis  equation  to  find  local  strains 
deltaeps(i)  =  deltasig(i)  /  E  +  2  *  (deltasig(i)  /  (2  *  xKf )  ) 
(1  /  xnf) 

eps(i)  =  lasteps  +  ((-1)  *  (flagl  +  i) )  *  deltaeps(i) 
lasteps  =  eps(i) 


END  IF 


'  pick  between  various  strain-life  eqns 
SELECT  CASE  menu2 . equation 

CASE  1  '  Morrow's  equation 

CALL  EQUATIONS 2 

CASE  2  '  Smith  Watson  Topper 

CALL  EQUATIONS 3 

CASE  3  ’  Manson  Halford 

CALL  EQUATIONS 4 
CASE  ELSE 
END  SELECT 


'  print  running  output  to  file 
IF  menu2 . option4  =  2  THEN  DATADUMP 
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IF  usedlife  >  1  THEN 
lasti  =  i 
EXIT  FOR 
END  IF 


NEXT  i 

IF  flag5  =  "OFF"  THEN  BEEP 

LOOP  UNTIL  usedlife  >=  1  OR  (menu2 . option4  =  2  AND  block  =  l) 
CLOSE  #7 
END  SUB 


D .  SUBROUTINE  DATADOMP 


REM  $ STATIC 

'  subroutine  to  write  output  data 

SUB  DATADUMP 

SHARED  menu2  AS  Options2 

IF  menu2 . inputs  =  "Stress"  THEN 

PRINT  #7,  USING  "  ###  ###  ###  ###.##  ####.## 

####.##  #.######  #.######  #.## - ";  i; 

stress (i);  deltaStress  (i) ;  deltasig(i);  sig(i); 
sigO(i);  deltaeps(i);  eps(i);  NNf(i) 

ELSE 

PRINT  #7,  USING  "  ###  #.###  #.###  ###.###  ###.##  ####.## 

####.##  #.######  #.######  i; 

strain (i ) ;  deltastrain (i) ;  deltaStress (i) ; 
deltasig(i);  sig(i) ;  sigO(i);  deltaeps(i); 
eps (i) ;  NNf (i) 

END  IF 

END  SUB 
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E .  SUBROUTINE  EQUATIONS! 


==========================  EQUATIONS1  =================== 

subroutine  to  find  delta  sigma  using  Newton's  method 
(based  on  Neuber's  Relation) 


SUB  EQUATIONSl 

DIM  Yprime  AS  DOUBLE 

Jeltasig(i)  =  1 
loopcount  =  0 


loopcount  =  loopcount  +  1 

Y  =  deltasig(i)  *  (deltasig(i)  /  (2  *  E)  +  (deltasig(i)  /  (2  *  xKf ) ) 

(1  /  xnf) )  -  ( Ktf  2)  *  deltaStress (i)  *  (deltaStress ( i )  /  (2  *  E) 

+  (deltaStress (i)  /  (2  *  xKf ) )  (1  /  xnf)) 

IF  ABS(Y)  >  .0000001  THEN 

Yprime  =  deltasig(i)  /  E  +  ({xnf  +  1)  /  xnf)  *  (deltasig(i)  /  (2  *  xKf)  ) 
(1  /  xnf) 

deltasig(i)  =  deltasig(i)  -  Y  /  Yprime 
END  IF 

LOOP  UNTIL  (ABS(Y)  <=  .0000001)  OR  (loopcount  =  10000) 

END  SUB 


F .  SUBROUTINE  EQUATIONS2 


'  subroutine  to  evaluate  Morrow's  Strain-Life  equation 


SUB  EQUATIONS 2 
REM  5DYNAMIC 

DIM  Yprime  AS  DOUBLE,  Y  AS  DOUBLE,  loopcount  AS  LONG 

NNf(i)  =  10000 
loopcount  =  0 


loopcount  =  loopcount  +  1 

Y  =  -deltaeps(i)  /  2  +  ( (sigff  -  sig0(i))  /  E)  *  (2  *  NNf(i))  b  +  epsff  * 
(2  *  NNf (i) )  *  c 

IF  ABS(Y)  >  .0000000001#  THEN 

Yprime  =  (b  *  (sigff  -  sig0(i))  /  E)  *  (2  A  b)  *  (NNf (i) )  *  (b  -  1)  +  c 
*  epsff  *  (2  *  c)  *  (NNf(i))  A  (c  -  1) 

IF  (Y  /  Yprime  <  NNf(i))  THEN 

NNf ( i)  =  NNf (i)  -  Y  /  Yprime 

ELSE 

NNf(i)  =  NNf (i)  /  2 

END  IF 
END  IF 
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LOOP  UNTIL  (ABS(Y)  <=.0000000001#)  OR  (NNf(i)  >  100000000)  OR  ( loopcount=10000 ) 
usedlife  =  usedlife  +  1  /  NNf(i) 


IF  loopcount  =  10000  THEN  PRINT  "*"  'indication  of  convergence  dificulties 

END  SUB 


G.  SUBROUTINE  EQUATIONS3 


REM  $ STATIC 

'  subroutine  to  evaluate  Smith  Watson  Topper  Strain-Life  equation 

SUB  EQUATIONS3 
REM  $ DYNAMIC 

DIM  Yprime  AS  DOUBLE 

NNf(i)  =  100 
loopcount  =  0 

DO 

loopcount  =  loopcount  +  l 

Y  =  -(sig0(i)  +  deltasig(i)  /  2)  *  (deltaeps(i)  /  2)  +  ( (sigff  A  2)  /  E)  *  (2 
*  NNf(i))  A  (2  *  b)  +  sigff  *  epsff  *  (2  *  NNf(i))  A  (b  +  c) 

IF  ABS(Y)  >  .0000000001#  THEN 

Yprime  =  (2  *  b  *  (sigff  '2)  /  E)  *  (2  A  (2  *  b) )  *  NNf(i)  A  (2  *  b  -  1) 
+  (b  +  c)  *  sigff  *  epsff  *  (2  A  (b  +  c) )  *  NNf(i)  A  (b  +  c  -  1) 

IF  (Y  /  Yprime  <  NNf(i))  THEN 

NNf(i)  =  NNf(i)  -  Y  /  Yprime  *  .5 
ELSE 

NNf  (i)  =  NNf  (i)  /  2 .5 
END  IF 
END  IF 

LOOP  UNTIL  (ABS(Y)  <=  . 0000000001#)  OR  (NNf(i)  >  100000000)  OR  (loopcount  = 
10000) 

usedlife  =  usedlife  +  1  /  NNf (i) 

IF  loopcount  =  10000  THEN  PRINT  'indication  of  convergence  dificulties 

END  SUB 
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H.  SUBROUTINE  EQUATIONS4 


'============================  EQUATIONS 4  ==================== 

'  subroutine  to  evaluate  Manson-Halford' s  Strain-Life  equation 

SUB  EQUATI0NS4 
REM  $ DYNAMIC 

DIM  Yprime  AS  DOUBLE 

NNf  (i)  =  100 
loopcount  =  0 


loopcount  =  loopcount  +  1 

Y  =  -deltaeps(i)  /  2  +  ( (sigff  -  sigO(i))  /  E)  *  (2  *  NNf(i))  b  +  epsff  * 
( (sigff  -  sigO(i))  /  sigff)  (b  /  c)  *  (2  *  NNf(i))  *  c 
IF  ABS(Y)  >  .0000000001#  THEN 

Yprime  =  (b  *  (sigff  -  sig0(i))  /  E)  *  (2  *  b)  *  NNf(i)  *  (b  -  1)  +  c  * 
epsff  *  ((sigff  -  sig0(i))  /  sigff)  (b  /  c)  *  (2  A  c)  *  NNf(i) 
A  (c  -  1) 

IF  (Y  /  Yprime  <  NNf(i))  THEN 

NNf(i)  =  NNf (i)  -  Y  /  Yprime  *  .5 
ELSE 

NNf ( i )  =  NNf ( i )  /  2 . 5 

END  IF 
END  IF 

LOOP  UNTIL  (ABS(Y)  <=.0000000001#)  OR  (NNf(i)  >  100000000)  OR  ( loopcount=10000 ) 
usedlife  =  usedlife  +  1  /  NNf(i) 

IF  loopcount  =  10000  THEN  PRINT  'indication  of  convergence  dificulties 

END  SUB 


I .  SUBROUTINE  EQUATIONS5 


REM  OSTATIC 

'  subroutine  to  find  the  far- field  stress  associated  with  the 
'  far-field  strain  inputs 

»  _  _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _  _ _ _ _ 

SUB  EQUATIONS 5 
DIM  Yprime  AS  DOUBLE 

deltaStress (i)  =  0 
loopcount  =  0 


IF  strain (i)  <  strainmax  THEN 
DO 

loopcount  =  loopcount  +  1 

Y  =  -deltastrain (i)  +  deltaStress (i)  /  E  +  2  *  (deltaStress  /  (2  *  xKf ) ) 

(1  /  xnf ) 
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IF  ABS (Y)  >  .0001  THEN 

Yprime  =  1/  E+  (1/  (2  *  xKf)  )  A  (1  /  xnf )  *  (2  /  xnf )  *  deltaStress 
(1  /  xnf  -  1) 

deltaStress (i)  =  deltaStress (i )  -  Y  /  Yprime 
END  IF 

LOOP  UNTIL  (ABS (Y)  <=  .0001)  OR  (loopcount  =  1000) 

'  calcalate  far  field  stress  from  change  in  stress 
stress(i)  =  laststress  +  ((-1)  *  (flagl  +  i) )  *  deltaStress (i) 
laststress  =  stress (i) 

ELSE 

DO 

loopcount  =  loopcount  +  l 

Y  =  -strain(i)  +  stress <i)  /  E  +  (stress  /  xKf)  (1  /  xnf) 

IF  ABS (Y)  >  .0001  THEN 

Yprime  =  1  /  E  +  (1  /  (xKf) )  *  (1  /  xnf)  *  (1  /  xnf)  *  stress  A  (1 
/  xnf  -  1) 

stress (i)  =  stress (i)  -  Y  /  Yprime 
END  IF 

LOOP  UNTIL  (ABS (Y)  <=  .0001)  OR  (loopcount  =  1000) 

'  cal.  farfield  stress  change  and  reset  strainmax 

strainmax  =  strain (i) 

deltaStress (i)  =  ABS (stress (i)  -  laststress) 
laststress  =  stress (i) 

END  IF 


END  SUB 


J.  SUBROUTINE  GetConfig 


DEFSNG  I-P 


Get  the  starting  number  of  lines  and  the  video  adapter. 


SUB  GetConfig  STATIC 

SHARED  InitRows  AS  INTEGER,  BestMode  AS  INTEGER,  Available  AS  STRING 

'  Assume  50  line  display  and  fall  through  error 

'  until  we  get  the  actual  number 

InitRows  =  50 

ON  ERROR  GOTO  RowErr 

LOCATE  InitRows,  1 

'  Assume  best  possible  screen  mode 
BestMode  =  VGA 
Available  =  "12789BCD" 

ON  ERROR  GOTO  VideoErr 

'  Fall  through  error  trap  until  a  mode  works 
SCREEN  BestMode 

'  If  EGA,  then  check  pages  to  see  whether  more  than  64K 
ON  ERROR  GOTO  EGAErr 

IF  BestMode  =  EGA2SG  THEN  SCREEN  8,  ,  1 
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ON  ERROR  GOTO  0 


'  Reset  text  mode 
SCREEN  0,  ,  0 

WIDTH  80,  25 

END  SUB 


K.  SUBROUTINE  HEADER 


DEFINT  I-P 


===========================  HEADER  =============================== 

subroutine  to  an  appropriate  output  header 


SUB  HEADER 


SHARED  menu2  AS  Options2 


IF  menu2. inputs 
PRINT  #7, 

ELSE 

PRINT  #7 ,  " 

END  IF 


=  "Stress"  THEN 

"index  Stress  deltaStress  deltasig 
deltaeps  eps  NNf  " 


index  strain  dstrn  deltaStress  deltasig 
deltaeps  eps  NNf  " 


sig 


sig 


END  SUB 


L .  SUBROUTINE  Klaxon 


DEFINT  H 

'  subroutine  activates  a  two  tone  alarm  until  a  key  is  pressed 


SUB  Klaxon  (Hi% ,  low%)  STATIC 

DO  WHILE  INKEY$  =  "" 

SOUND  Hi,  5 
SOUND  low,  5 
LOOP 

END  SUB 


sigO 

sigO 
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M.  SUBROUTINE  LOADER 


REM  $ DYNAMIC 
DEFSNG  H-P 

' ==============================  LOADER  =============================== 

'  subroutine  to  read  in  the  array  of  loads  from  the  file  "STRESS.DAT" 
'  or  "STRAIN.DAT"  as  appropriate 

'  also  --  "top"  is  set  to  the  highest  index  in  the  load  array  and 

'  --  "flagl"  is  set  to  indicate  if  the  loads  are  initially 

'  increasing  or  decreasing 


SUB  LOADER 

SHARED  menu2  AS  0ptions2 
OPEN  inputfile  FOR  INPUT  AS  #10 

j  =  1 

IF  menu2 . inputs  =  "Stress"  THEN 

DO  UNTIL  EOF (10) 

INPUT  #10,  stress  (j) 
3=1+1 
LOOP 


top  =  j  -  1 


IF  stress  (1) 
flagl  = 

ELSE 

flagl  = 

END  IF 

ELSE 

DO  UNTIL  EOF (10) 

INPUT  #10,  strain (j) 

j=j+l 

LOOP 

'  set  "top"  to  the  number  of  loads: 

top  =  j  -  l 

'  set  "flagl":  (same  as  above) 

IF  strain (1)  >  0  THEN 
flagl  =  1 
ELSE 

flagl  =  0 
END  IF 

END  IF 

CLOSE  #10 

END  SUB 


set  "top"  to  the  number  of  loads: 


set  "flagl":  (indicating  first  cycle 
>  0  THEN  '  in  tension  or  compression) 

1 

0 
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N.  SUBROUTINE  Loa.dma.terial 


REM  $STATIC 
DEFINT  I-P 

' ===========================  Loadmatenal  =================== 

'  finds  the  selected  material  in  the  data  base  and  loads  the 
'  material  values  for  the  selected  material 

'  also  calculates  any  critical  missing  values  it  can 


SUB  Loadmaterial  (index) 

'OPEN  "A:\MAT.DAT"  FOR  INPUT  AS  #3 
'OPEN  "B:\MAT.DAT"  FOR  INPUT  AS  #3 
OPEN  "MAT. DAT"  FOR  INPUT  AS  #3 

INPUT  #3,  count 

DO  UNTIL  index  =  matindex 

INPUT  #3,  matindex,  matname$,  Su,  Sy,  Syf,  K,  Kf,  n,  nf,  epsf,  epsff,  sigf, 
sigff,  b,  c,  Sf,  SfSu,  E 
LOOP 

CLOSE  #3 

IF  n  =  0  AND  b  <>  0  AND  c  <>  0  THEN  n  =  b  /  c 

IF  nf  =  0  AND  b  <>  0  AND  c  <>  0  THEN  nf  =  b  /  c 

IF  K  =  0  AND  n  <>  0  AND  epsf  <>  0  THEN  K  =  sigf  /  (epsf  A  n) 

IF  Kf  =  0  AND  nf  <>  0  AND  epsff  o  0  THEN  Kf  =  sigff  /  (epsff  A  nf) 

END  SUB 


0.  FUNCTION  LOGIO 


'  ============================  LOGIO  ================== 

'  Returns  the  base  10  logarithm  for  the  specified  value 

FUNCTION  LOGIO i  (value) 

LOGIO  =  LOG (value)  /  LOG (10) 

END  FUNCTION 
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P .  SUBROUTINE  MATMENU 


========================  MATMENU  =============================== 

subroutine  to  update  the  screen  when  entering  or  viewing  materials 
in  the  materials  database 


SUB  MATMENU 


CONST  ENTER  =  13,  ESCAPE  =27 

CONST  DOWNARROW  =  80,  UPARROW  =  72,  LEFTARROW  =  75,  RIGHTARROW  =77 
CONST  COL1  =  10,  COL2  =  50,  ROW  =  7 
CONST  COL3  =  7,  COL4  =  62 
CONST  Fields  =  14 

SHARED  Menu  AS  Opt ions 1 

DIM  forceunit  AS  STRING,  Fid  AS  INTEGER 
CLS 

IF  f lag3  =  0  THEN 

'  Display  key  instructions 
LOCATE  1,  COL1 

PRINT  "UP  .  Move  to  next  field" 

LOCATE  2,  COL1 

PRINT  "DOWN  .  Move  to  previous  field" 

LOCATE  3,  COL1 

PRINT  "LEFT  or  RIGHT  ....  Enter  a  new  value" 

LOCATE  4,  COL1 

PRINT  "ENTER  . . .  Return  with  current  values" 

LOCATE  5,  COL1 

PRINT  "ESCAPE  ....  Quit  material  data  entry" 

ELSE  'display  material's  name  at  top  of  screen 
LOCATE  2,  COL1 

PRINT  "ESCAPE  ....to  return  to  previous  previous  menu" 

LOCATE  5,  COL1 

PRINT  USING  "Data  base  parameters  for:  ;  matname 

END  IF 

IF  (Menu .material  <=  51)  THEN 
forceunit  =  "ksi) " 

ELSE 

forceunit  =  "MPa) " 

END  IF 


'  Display  fields 

LOCATE  ROW,  COL3 :  PRINT  "Ultimate  Strength 
LOCATE  ROW,  COL4 :  PRINT  USING  "[  ####  ]";  Su; 

LOCATE  ROW  +  1,  COL3 :  PRINT  "Yield  Strength 

LOCATE  ROW  +  1,  COL4 :  PRINT  USING  "[  ####  ] " ;  Sy; 

LOCATE  ROW  +  2,  COL3 :  PRINT  "Fatigue  Yield  Strength 

LOCATE  ROW  +  2,  COL4 :  PRINT  USING  "[  ####  ] " ;  Syf ; 


(Su  in  " ; 
forceunit 


(Sy  in 
forceunit ; 


(Sy'  in 
forceunit ; 
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LOCATE 

ROW 

+ 

3, 

COL3  : 

PRINT  "Strength  Coeficient 

(K  in  "  ; 
f orceunit ; 

LOCATE 

ROW 

+ 

3, 

COL4  : 

PRINT  USING  "[  ####  ]";  K; 

LOCATE 

ROW 

+ 

4, 

COL3  : 

PRINT  "Cyclic  Strength  Coefficient 

(K'  in  "  ; 
forceunit ; 

LOCATE 

ROW 

+ 

4  , 

COL4  : 

PRINT  USING  " [  ] " ;  Kf; 

LOCATE 

ROW 

+ 

5, 

COL3  : 

PRINT  "Strain  Hardening  Exponent 

(n)  "  ; 

LOCATE 

ROW 

+ 

5, 

COL4  : 

PRINT  USING  "[  #.##  ]";  n; 

LOCATE 

ROW 

+ 

6, 

COL3  : 

PRINT  "Cyclic  Strain  Hardening  Exponent 

(n' ) "  ; 

LOCATE 

ROW 

+ 

6, 

COL4  : 

PRINT  USING  "[#.##]";  nf; 

LOCATE 

ROW 

+ 

7, 

COL3  : 

PRINT  "Ductility  Coefficient 

(epsilon  f )  " 

LOCATE 

ROW 

+ 

7, 

COL4  : 

PRINT  USING  "[  ]";  epsf; 

LOCATE 

ROW 

+ 

8 ,  COL3  : 

PRINT  "Fatigue  Ductility  Coefficient 

(epsilon  f '  )  " 

LOCATE 

ROW 

+ 

8, 

COL4  : 

PRINT  USING  "[  epsff; 

LOCATE 

ROW 

+ 

9, 

COL3  : 

PRINT  "Strength  Coefficient 

(sigma  f  in  " 
forceunit 

LOCATE 

ROW 

+ 

9, 

COL4  : 

PRINT  USING  "[  ####  ]";  sigf; 

LOCATE 

ROW 

+ 

o 

H 

COL3  : 

PRINT  "Fatigue  Strength  Coefficient  ( 

sigma  f'  in  " 
forceunit 

LOCATE 

ROW 

+ 

10, 

COL4 

:  PRINT  USING  ”[  ####  ]";  sigff; 

LOCATE 

ROW 

+ 

ii. 

COL3 

:  PRINT  "Fatigue  strength  Exponent 

(b)  " : 

LOCATE 

ROW 

+ 

ii, 

COL4 

-  1:  PRINT  USING  "[  ##.##  ] " ;  b; 

LOCATE 

ROW 

+ 

12, 

COL3 

:  PRINT  "Fatigue  Ductility  Exponent 

(c)  "  ; 

LOCATE 

ROW 

+ 

12, 

COL4 

-  1:  PRINT  USING  "[  ##.##  ]";  c; 

LOCATE 

ROW 

+ 

13, 

COL3 

:  PRINT  "Endurance  Strength 

(Sf  in  " ; 
forceunit 

LOCATE 

ROW 

13, 

COL4 

-  1:  PRINT  USING  "[  #####  Sf; 

LOCATE 

ROW 

+ 

14, 

COL3 

:  PRINT  "Modulus  of  Elasticity 

(E  in  "  ; 
forceunit 

LOCATE 

ROW 

+ 

14, 

COL4 

-  1:  PRINT  USING  "[  #####  ]  " ;  E ; 

END  SUB 
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Q.  SUBROUTINE  NeuberKf 


'  subroutine  to  calculate  the  fatigue  stress  concentration  factor 
'  based  on  Neuber's  method 


SUB  NeuberKf 

CALL  TYPINCEnter  notch  sensitivity  factor  (q;  0  if  unknown)",  q) 

IF  q  <>  0  THEN 

Ktf  =  1  +  q  *  (Kt  -  1) 

ELSE 

CALL  TYPINCEnter  the  notch  root  radius  (r)",  r) 

CALL  TYPINCEnter  the  appropriate  material  constant  (rho)",  rho) 
Ktf  =  1  +  (Kt  -  1)  /  (1  +  (rho  /  r)  A  .5) 

END  IF 

END  SUB 


R.  SUBROUTINE  NEWMAT 


========================  NEWMAT  =========================== 

subroutine  to  manually  enter  material  constants  and 

enter  new  materials  into  the  materials  database 


SUB  NEWMAT  (Ky) 

CONST  ENTER  =  13,  ESCAPE  =  27 

CONST  DOWNARROW  =  80,  UPARROW  =  72,  LEFTARROW  =  75,  RIGHTARROW  =  77 
CONST  COL1  =  10,  C0L2  =  50,  ROW  =  7 
CONST  COL3  =  7,  COL4  =62 
CONST  Fields  =  14 


'  conversion  constant  from  "ksi"  to  "MPa" 
CONST  ab  =  6.89476 


SHARED  Menu  AS  Optionsl,  save$ 

DIM  forceunit  AS  STRING,  Fid  AS  INTEGER 


IF  (Menu .material  =  1)  THEN 
forceunit  =  "ksi) " 

ELSE 

forceunit  =  "MPa) " 

END  IF 

CALL  MATMENU 


DO 


Update  field  values  and  position  based  on  keystrokes 

'  Put  cursor  on  field 
LOCATE  ROW  +  Fid,  COL4  +  2 
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DO 


Get  a  key  and  strip  null  off  if  it's  an  extended  code 


Key$  =  INKEY$ 

LOOP  WHILE  Key$  =  " " 

Ky  =  ASC (RIGHTS (Key$,  1) ) 

SELECT  CASE  Ky 
CASE  ENTER 

'  Chech  for  all  appropriate  data  parameters 
'  before  allowing  return  to  main  option  menu 
IF  b  =  0  OR  c  =  0  OR  sigff  =  0  OR  epsff  THEN 

LOCATE  23,  10:  PRINT  "Insufficient  material  parameters  were  entered" 
END  IF 

CASE  UP ARROW,  DOWNARROW 
'  Adjust  field  location 

IF  Ky  =  DOWNARROW  THEN  Inc  =  1  ELSE  Inc  =  -1 
Fid  -  Rotated(0,  Fields,  Fid,  Inc) 

CASE  RIGHTARROW,  LEFTARROW 
'  Adjust  field  value 

IF  Ky  =  RIGHTARROW  THEN  Inc  =  1  ELSE  Inc  =  -1 
SELECT  CASE  Fid 

CASE  0 

'  Ultimate  Strength 

CALL  TYPIN ( "Enter  Ultimate  Strength  (Su)",  Su) 

CALL  MATMENU 

CASE  1 

'  Yield  Strength 

CALL  TYPIN ("Enter  Yield  Strength  (Sy)",  Sy) 

CALL  MATMENU 

CASE  2 

'  Fatigue  Yield  Strength 

CALL  TYPIN ( "Enter  Fatigue  Yield  Strength  (Sy')",  Syf) 

CALL  MATMENU 

CASE  3 

'  Strength  Coeficient 

CALL  TYPIN ( "Enter  Strength  Coeficient  (K)",  K) 

CALL  MATMENU 

CASE  4 

'  Cyclic  Strength  Coeficient 

CALL  TYPIN ("Enter  Cyclic  Strength  Coeficient  (K')",  Kf) 

CALL  MATMENU 

CASE  5 

'  Cyclic  Strength  Coeficient 

CALL  TYPIN( "Strian  Hardening  Exponent  (n) " ,  n) 

CALL  MATMENU 

CASE  6 

'  Cyclic  Strian  Hardening  Exponent 

CALL  TYPIN( "Cyclic  Strian  Hardening  Exponent  (n')",  nf) 
CALL  MATMENU 
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CASE  7 

'  Ductility  Coefficient  (epsilon  f) 

CALL  TYPIN ( "Ductility  Coefficient  (epsilon  f)",  epsf) 

CALL  MATMENU 

CASE  8 

'  Fatigue  Ductility  Coefficient  (epsilon  f') 

CALL  TYPIN ( "Fatigue  Ductility  Coefficient  (epsilon  f')", 
epsf f ) 

CALL  MATMENU 
CASE  9 

'  Strength  Coefficient  (sigma  f) 

CALL  TYPIN ( "Strength  Coefficient  (sigma  f)",  sigf) 

CALL  MATMENU 

CASE  10 

'  Fatigue  Strength  Coefficient  (sigma  f') 

CALL  TYPIN ("Fatigue  Strength  Coefficient  (sigma  f')",  sigff) 
CALL  MATMENU 

CASE  11 

'  Fatigue  Strength  Exponent  (b) 

CALL  TYPIN ("Fatigue  Strength  Exponent  (b) " ,  b) 

CALL  MATMENU 

CASE  12 

'  Fatigue  Ductility  Exponent  (c) 

CALL  TYPIN ( "Fatigue  Ductility  Exponent  (c) " ,  c) 

CALL  MATMENU 

CASE  13 

'  Endurance  Strength  (Sf) 

CALL  TYPIN ("Endurance  Strength  (Sf)",  Sf) 

CALL  MATMENU 

CASE  14 

'  Modulus  of  Elasticity  (E) 

CALL  TYPIN ("Modulus  of  Elasticity  (E)”,  E) 

CALL  MATMENU 

CASE  ELSE 
END  SELECT 

CASE  ELSE 
END  SELECT 

'  exit  do  loop  and  continue  with  subroutine  if  ENTER 
LOOP  UNTIL  (Ky  =  ENTER  OR  Ky  =  ESC) 


CLS 

LOCATE  13,  5:  PRINT  "Enter  Material's  Name  (up  to  30  charactors) :  " 
LOCATE  13,  52:  INPUT  matname$ 

LOCATE  15,  5:  PRINT  "Save  this  material  in  material  data  base  (Y/N) :  " 
LOCATE  15,  54:  INPUT  save$ 
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'  progam  segment  to  save  new  materials  in  the  material  data  base 
IF  saves  =  "Y"  OR  saveS  =  "y"  THEN 

'OPEN  "A:\MAT.DAT"  FOR  APPEND  AS  #3 
'OPEN  "B:\MAT.DAT"  FOR  APPEND  AS  #3 
OPEN  "MAT. DAT"  FOR  APPEND  AS  #3 

matcount  =  matcount  +  1 
IF  forceunit  =  "ksi) "  THEN 
matindex  =  matcount 

WRITE  #3,  matindex,  matname$,  Su,  Sy,  Syf,  K,  Kf,  n,  nf,  epsf,  epsff, 
sigf,  sigff,  b,  c,  Sf,  SfSu,  E 

WRITE  S3,  matindex  +  50,  matname$,  Su  *  ab,  Sy  *  ab,  Syf  *  ab,  K  *  ab, 

Kf  *  ab,  n,  nf,  epsf,  epsff,  sigf  *  ab,  sigff  *  ab,  b,  c,  Sf 

*  ab,  SfSu,  E  *  ab 

ELSE 

matindex  =  matcount  +  50 

WRITE  #3,  matindex,  matname$,  Su,  Sy,  Syf,  K,  Kf,  n,  nf,  epsf,  epsff, 
sigf,  sigff,  b,  c,  Sf,  SfSu,  E 

WRITE  S3,  matindex  -  50,  matname$,  Su  /  ab,  Sy  /  ab,  Syf  /  ab,  K  /  ab, 

Kf  /  ab,  n,  nf,  epsf,  epsff,  sigf  /  ab,  sigff  /  ab,  b,  c,  Sf 

/  ab,  SfSu,  E  /  ab 

END  IF 
CLOSE  S3 

'OPEN  "A:\NEWCOUNT.DAT"  FOR  OUTPUT  AS  #4 
'OPEN  "B:\NEWCOUNT.DAT"  FOR  OUTPUT  AS  S4 
OPEN  " NEWCOUNT.DAT"  FOR  OUTPUT  AS  #4 
newstuff  =  newstuff  +  1 
WRITE  #4 ,  newstuff 
CLOSE  #4 

END  IF 

CLS 

CALL  UPDATEMENU 
Ky  =  0 

END  SUB 


S .  SUBROUTINE  OPTMENU1 

REM  $ DYNAMIC 

'  ===========================  Option  Menu  1  ========================== 

'  Def ine/swtich  the  user  selectable  functions  for  the  material  and 

'  stress  concentration  factors 

SUB  OPTMENU1  STATIC 

SHARED  VC  AS  Config,  Menu  AS  Optionsl,  Available  AS  STRING,  Fields  AS  INTEGER 
SHARED  Fid  AS  INTEGER 

'  Constants  for  key  codes  and  column  positions 
CONST  ENTER  =  13,  ESCAPE  =27,  FI  =  59 

CONST  DOWNARROW  =  80,  UPARROW  =  72,  LEFTARROW  =  75,  RIGHTARROW  =77 
CONST  COL1  =  10,  COL2  =  50,  ROW  =  9 
CONST  COL3  =  7,  COL4  =42 
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'OPEN  "A:\MAT.DAT"  FOR  INPUT  AS  #3 
'OPEN  "B:\MAT.DAT"  FOR  INPUT  AS  #3 
OPEN  "MAT. DAT"  FOR  INPUT  AS  #3 
INPUT  #3,  matcount 
CLOSE  #3 

'OPEN  "A:\NEWCOUNT.DAT”  FOR  INPUT  AS  #4 
'OPEN  "B:\NEWCOUNT.DAT"  FOR  INPUT  AS  #4 
OPEN  "NEWCOUNT.DAT"  FOR  INPUT  AS  #4 
INPUT  #4,  newstuff 
CLOSE  #4 

matcount  =  matcount  +  newstuff 
CALL  Loadmaterial (Menu. material) 


'  Block  cursor 

LOCATE  ROW,  COL1 ,  1,  1,  12 

CALL  UPDATEMENU 

'  Skip  field  10  if  there's  only  one  value 

IF  LEN ( Avai lable$ )  =  l  THEN  Fields  =  8  ELSE  Fields  =  10 


DO 


Update  field  values  and  position  based  on  keystrokes 

'  Put  cursor  on  field 
LOCATE  ROW  +  Fid,  COL4  +  2 


'  Get  a  key  and  strip  null  off  if  it's  an  extended  code 
DO 

Key$  =  INKEY$ 

LOOP  WHILE  Key$  =  "" 

Ky  =  ASC (RIGHT$ (Key$ ,  1)) 


SELECT  CASE  Ky 
CASE  ESCAPE 

'  End  program 
CLS  :  END 

CASE  FI  '  review  material  parameters 

flag3  =  1 
CALL  MATMENU 
flag3  =  0 
DO 

DO 

Key2$  =  INKEY$ 

LOOP  WHILE  Key2$  =  "" 

Ky2  =  ASC (RIGHT$ (Key2$ ,  1)) 

LOOP  UNTIL  Ky2  =  ESCAPE 
CLS 

CALL  UPDATEMENU 
CASE  ENTER 

'  runs  material  parameter  entry  subroutine 
IF  Menu. material  =  1  OR  Menu .material  =  51  THEN  NEWMAT  (Ky) 
CASE  UPARROW,  DOWNARROW 
'  Adjust  field  location 

IF  Ky  =  DOWNARROW  THEN  Inc  =  2  ELSE  Inc  =  -2 
Fid  =  Rotated (0,  Fields,  Fid,  Inc) 
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CASE  RIGHTARROW,  LEFTARROW 
'  Adjust  field  value 

IF  Ky  =  RIGHTARROW  THEN  Inc  =  1  ELSE  Inc  =  -1 
SELECT  CASE  Fid 

CASE  0 

'  Units 

IF  Menu. units  =  "Brit"  THEN 
Menu. units  =  "  SI  " 

Menu. material  =  Menu. material  +  50 
CALL  Loadmaterial (Menu. material) 

CALL  UPDATEMENU 
ELSE 

Menu. units  =  "Brit" 

Menu .material  =  Menu. material  -  50 
CALL  Loadmaterial (Menu. material) 

CALL  UPDATEMENU 
END  IF 

CASE  2 

'  Material  Selection 
IF  Menu. units  =  "Brit"  THEN 

Menu .material  =  Rotated (1,  matcount,  Menu. material ,  Inc) 
CALL  Loadmaterial (Menu .material ) 

PRINT  USING  matname 

ELSE 

Menu. material  =  Rotated(51,  matcount  +  50,  Menu .material , 

Inc) 

CALL  Loadmaterial (Menu. material) 

PRINT  USING  "Sc";  matname 
END  IF 

CASE  4 

'  Stress  Consentration  Factor  Kt 

CALL  TYPIN( "Enter  stress  concentration  factor  (Kt)  :  ",  Kt) 

CALL  UPDATEMENU 

CASE  6 

'  Fatigue  Stress  Concentration  Calculation  Method 
Menu.Ktf  =  Rotated (1,  3,  Menu.Ktf,  Inc) 

SELECT  CASE  Menu.Ktf 
CASE  1 

'  manual  entry' 

option4  =  "Manually  Entered  " 

LOCATE  ROW  +  8,  COL4 ;  PRINT  USING 

LOCATE  ROW  +  6,  C0L4  +  2:  PRINT  USING  option4 

CASE  2 

'  Neubers  Method 

option4  =  "Neuber's  Method  " 

LOCATE  ROW  +  8,  COL4  :  PRINT  USING  "  £. "  ;  "[---] 

LOCATE  ROW  +  6,  COL4  +  2:  PRINT  USING  "Sc";  option4 
CASE  3 

'  Peterson  Method 

option4  =  "Peterson's  Method" 

LOCATE  ROW  +  8,  COL4 :  PRINT  USING 

LOCATE  ROW  +  6,  C0L4  +  2:  PRINT  USING  "Sc";  option4 
CASE  ELSE 
END  SELECT 
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CASE  8 

'  Fatigue  Stress  Concentration  Factor  displayed 
SELECT  CASE  Menu.Ktf 
CASE  1 

'  manual  entry' 

CALL  TYPIN ( "Enter  fatigue  stress  concentration  factor 
(Ktf) :  ",  Ktf ) 

CALL  UPDATEMENU 
CASE  2 

'  Neubers  Method 
CALL  NeuberKf 
CALL  UPDATEMENU 
CASE  3 

'  Peterson  Method 
CALL  PetersonKf 
CALL  UPDATEMENU 
CASE  ELSE 
END  SELECT 

CASE  10 

'  Available  screen  modes 
i  =  INSTR  (Available$,  HEX$  (VC  .  Scm)  ) 
i  =  Rotated(l,  LEN (Available$) ,  i,  Inc) 

VC. Sera  =  VAL ( " &h "  +  MID$ (Available? ,  i,  1) ) 

PRINT  USING  VC.Scrn 

CASE  ELSE 
END  SELECT 

CASE  ELSE 
END  SELECT 


'  if  cluase  to  ensure  selection  of  a  fatigue  stress  concentation  factor 
IF  Ky  =  ENTER  AND  Ktf  =  0  THEN 
Fid  =  8 

CALL  TYPINCEnter  fatigue  stress  concentration  factor  (Ktf):  ",  Ktf) 
CALL  UPDATEMENU 
Ky  =  0 
END  IF 


'  Return  to  main  program  if  ENTER 

LOOP  UNTIL  (Ky  =  ENTER  AND  Ktf  <>  0  AND  matname  <>  "new  material") 


END  SUB 
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T .  SUBROUTINE  0PTMENU2 


=========================  Option  Menu  2  ======================== 

Def ine/swtich  the  user  selectable  functions  for  processing  options 


SUB  0PTMENU2  STATIC 

SHARED  VC  AS  Config,  menu2  AS  0ptions2,  Available  AS  STRING,  Fields  AS  INTEGER 
SHARED  Fid  AS  INTEGER 

'  Constants  for  key  codes  and  column  positions 
CONST  ENTER  =  13,  ESCAPE  =27,  F2  =  60 

CONST  DOWNARROW  =  80,  UPARROW  =  72,  LEFTARROW  =  75,  RIGHTARROW  =77 
CONST  COL1  =  10,  COL2  =  50,  ROW  =  9 
CONST  COL3  =  7,  COL4  =42 

'  Return  cursor  to  menu  top 
Fid  =  0 

CALL  UPDATEMENU2 

'  Skip  field  10  if  there's  only  one  value 

IF  LEN (Available$)  =  1  THEN  Fields  =  8  ELSE  Fields  =  10 


DO 


Update  field  values  and  position  based  on  keystrokes 

'  Put  cursor  on  field 
LOCATE  ROW  +  Fid,  COL4  +  2 


'  Get  a  key  and  strip  null  off  if  it's  an  extended  code 
DO 

Key$  =  INKEY$ 

LOOP  WHILE  Key$  =  "" 

Ky  =  ASC(RIGHT$ (Key$,  1)) 

SELECT  CASE  Ky 
CASE  ESCAPE 

'  End  program 
CLS  :  END 
CLS  :  f lag2  =  0 

CASE  F2  '  changes  sound  to  condition  show  in  flags 

IF  flags  =  "ON"  THEN  '  (sound  is  opposite  of  flag) 
f lags  =  "OFF" 

ELSE 

flags  =  "ON" 

END  IF 
UPDATEMENU2 

CASE  UPARROW,  DOWNARROW 
'  Adjust  field  location 

IF  Ky  =  DOWNARROW  THEN  Inc  =  2  ELSE  Inc  =  -2 
Fid  =  Rotated(0,  Fields,  Fid,  Inc) 

CASE  RIGHTARROW,  LEFTARROW 
'  Adjust  field  value 

IF  Ky  =  RIGHTARROW  THEN  Inc  =  1  ELSE  Inc  =  -1 
SELECT  CASE  Fid 
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CASE  0 

'  type  of  inputs  (stress  vs  strain) 

IF  menu2 . inputs  =  "Stress"  THEN 
menu2. inputs  =  "Strain" 
inputfile  =  " STRAIN . DAT" 

ELSE 

menu2. inputs  =  "Stress" 
inputfile  =  " STRESS . DAT" 

END  IF 
UPDATEMENU2 

CASE  2 

'  Equation  Selection 

menu2 .equation  =  Rotatedd,  3,  menu2 . equation,  Inc) 
SELECT  CASE  menu2 . equation 
CASE  1 

'  Morrow's  equation 
equationname  =  "Morrow's  equation" 

CASE  2 

'  Smith  Watson  Topper 
equationname  =  " Smith- Watson-Topper" 

CASE  3 

'  Manson  Halford 
equationname  =  "Manson- Hal ford" 

CASE  ELSE 
END  SELECT 
UPDATEMENU2 

CASE  4 

'  variable /fixed  n'  and  K' 

menu2  ,option3  =  Rotatedd,  4,  menu2  ,option3 ,  Inc) 
SELECT  CASE  menu2 . option3 
CASE  1 

'  fixed  n  and  K 

nKtype  =  "  Fixed  n'  and  K'  " 

CASE  2 

'  variable  n'  and  K' 

nKtype  =  "Variable  n'  and  K' " 

CASE  3 

’  variable  n'  and  fixed  K' 

nKtype  =  "Variable  n'  and  fixed  K'  " 

CASE  4 

'  fixed  n'  and  variable  K' 
nKtype  =  "Fixed  n'  and  variable  K' " 

CASE  5 

'  Experimental  -  not  used 
nKtype  =  "Experimental" 

CASE  ELSE 
END  SELECT 
UPDATEMENU2 

CASE  6 

'  calculation  type 

menu2  . option4  =  Rotatedd,  3,  menu2  . option4 ,  Inc) 
SELECT  CASE  menu2 . option4 
CASE  1 

'  blocks  to  failure 

calctype  =  "Load  blocks  to  failure" 

CASE  2 

'  single  block  effects 
calctype  =  "Single  block  effects" 
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CASE  3 

'  batch  process 
calctype  =  "Batch  Process" 
flag4  =  l 
CASE  ELSE 
END  SELECT 
UPDATEMENU2 

CASE , 8 

'  input  data  file  name 

CALL  TYPINSTRING (" Input  data  file's  name:",  inputfile) 
UPDATEMENU2 

CASE  10 

'  output  data  file  name 

CALL  TYPINSTRING ( "Output  data  file's  name:",  outputfile) 
UPDATEMENU2 

CASE  ELSE 
END  SELECT 

CASE  ELSE 
END  SELECT 


'  if  cluase  to  ensure  selection  of  a  fatigue  stress  concentation  factor 
IF  Ky  =  ENTER  AND  Ktf  =  0  THEN 
Fid  =  8 

CALL  TYPIN( "Enter  fatigue  stress  concentration  factor  (Ktf):  ",  Ktf) 
CALL  UPDATEMENU 
Ky  =  0 
END  IF 


'  Return  to  main  program  if  ENTER 

LOOP  UNTIL  (Ky  =  ENTER  AND  Ktf  <>  0  AND  matname  <>  "new  material")  OR  (Ky  = 

ESCAPE) 


CLS 

END  SUB 
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U.  SUBROUTINE  OUTPUTER 


REM  $ STATIC 

'  subroutine  writes  to  an  output  file  compution  parameters 

'  and  results 


SOB  OUTPUTER 


OPEN  outputfile  FOR  APPEND  AS  #11 


WRITE  #11, 
WRITE  #11, 
WRITE  #11, 
WRITE  #11, 
WRITE  #11, 
WRITE  #11, 
WRITE  #11, 
WRITE  #11, 
WRITE  #11, 
WRITE  #11, 
WRITE  #11, 
WRITE  #11, 
WRITE  #11, 


equationname 

nKtype 

calctype 

"input  file:",  inputfile 
"output  file:",  outputfile 
"blocks:",  block 
"i  counter:",  lasti 
"reversal  count:",  NNf count 
"life  factor:",  usedlife 


CLOSE  #11 


END  SUB 


V.  SUBROUTINE  PetersonKf 


DEFSNG  I - P 

'  subroutine  to  calculate  the  fatigue  stress  concentration  factor 
'  based  on  Peterson's  method 


SUB  PetersonKf 

CALL  TYPIN( "Enter  notch  sensitivity  factor  (q;  0  if  unknown)",  q) 
IF  q  <>  0  THEN 

Ktf  =  1  +  q  *  (Kt  -  1) 

ELSE 

CALL  TYPIN( "Enter  the  notch  root  radius  (r)",  r) 

CALL  TYPIN( "Enter  the  appropriate  material  constant  (a)",  a) 
Ktf  =  1  +  (Kt  -  1)  /  (1  +  a  /  r) 

END  IF 
END  SUB 
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W.  FUNCTION  Rotated 


DEFINT  I-P 

'  Returns  the  present  value  adjusted  by  Inc  and  rotated  if  necessary 
'  so  that  it  falls  within  the  range  of  Lower  and  Upper. 

FUNCTION  Rotated  (Lower,  Upper  AS  INTEGER,  present,  Inc) 

'  Calculate  the  next  value 
present  =  present  +  Inc 

'  Handle  special  cases  of  rotating  off  top  or  bottom 
IF  present  >  Upper  THEN  present  =  Lower 
IF  present  <  Lower  THEN  present  =  Upper 
Rotated  =  present 

END  FUNCTION 


X.  SUBROUTINE  TYPIN 


' _===___======_===_=_==========  TYPIN  =========================== 

'  subroutine  to  update  the  value  of  a  specific  variable  through 

'  keyboard  entry  (real  number  variable) 

SUB  TYPIN  (cstring$ ,  valu) 

LOCATE  23,  10:  PRINT  cstring$ 

LOCATE  23,  58:  INPUT  valu 
LOCATE  23,  10:  PRINT  " 

END  SUB 


Y .  SUBROUTINE  TYPINSTRING 


'  subroutine  to  update  the  value  of  a  specific  string  variable  through 
'  keyboard  entry 

SUB  TYPINSTRING  (cstring$,  stringname$) 

LOCATE  23,  10:  PRINT  cstring$ 

LOCATE  23,  58:  INPUT  stringname$ 

LOCATE  23,  10:  PRINT  " 

END  SUB 
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Z .  SUBROUTINE  UPDATEMENO 


======================  UPDATEMENU  ============ 

subroutine  to  update  the  options  menu  (menu  1) 


SUB  UPDATEMENU 

SHARED  VC  AS  Config,  Menu  AS  Optionsl,  Available  AS  STRING,  Fields  AS  INTEGER, 
Fid  AS  INTEGER 

'  Constants  for  key  codes  and  column  Dositions 
CONST  C0L1  =  10,  COL2  =  50,  ROW  =  9 
CONST  COLS  =  7,  COL4  =42 


'  Display  key  instructions 
LOCATE  1,  COL1 

PRINT  "UP  .  Move  to  next  field" 

LOCATE  2,  COL1 

PRINT  "DOWN  .  Move  to  previous  field" 

LOCATE  3,  COL1 

PRINT  "LEFT/RIGHT  ....  Change  field  up/down" 
LOCATE  4,  COL1 

PRINT  "FI  ....  Display  matrial's  parameters" 
LOCATE  5,  COL1 

PRINT  "ENTER  ....  Start  with  current  values" 
LOCATE  6,  COL1 

PRINT  "ESCAPE  .  Quit  Program" 


'  Display  fields 

LOCATE  ROW,  COL3 :  PRINT  "Type  of  units  (SI  or  British)" 
LOCATE  ROW,  COL4 :  PRINT  USING  "[  &  ]";  Menu. units; 


LOCATE 

ROW 

+ 

2, 

COL3  : 

PRINT 

"Material" ; 

LOCATE 

ROW 

+ 

2, 

COL4  : 

PRINT 

USING  " [ 

&  ]  " ;  matname ; 

LOCATE 

ROW 

+ 

4, 

COL  3  : 

PRINT 

"Stress 

concentration  factor 

LOCATE 

ROW 

+ 

4, 

COL4  : 

PRINT 

USING  ” t 

##.###  ]";  Kt; 

LOCATE 

ROW 

+ 

6, 

COL3  : 

PRINT 

"Method 

to  calculate  Kf"; 

LOCATE 

ROW 

+ 

6, 

COL4  : 

PRINT 

USING  " [ 

Sc  ]  "  ;  option4 

LOCATE 

ROW 

+ 

8, 

COL3  : 

PRINT 

"Fatigue 

stress  cone,  factor 

LOCATE 

ROW 

+ 

8, 

COL4  : 

PRINT 

USING  " t 

##.###  ]";  Ktf; 

LOCATE 

ROW 

+ 

10, 

COL3 

:  PRINT  "Screen 

Mode  " ; 

LOCATE 

ROW 

+ 

10, 

COL4 

:  PRINT  USING  ” 

[  ##  ] " ;  VC . Scrn 

END  SUB 
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AA.  SUBROUTINE  OPDATEMENU2 


======================  UPDATEMENU2  =========== 

subroutine  to  update  the  options  menu  (menu  2) 


SUB  UPDATEMENU2 

SHARED  VC  AS  Config,  menu2  AS  0ptions2,  Available  AS  STRING,  Fields  AS  INTEGER, 
Fid  AS  INTEGER 

'  Constants  for  key  codes  and  column  positions 
CONST  CO LI  =  10,  COL2  =  50,  ROW  =  9 
CONST  COL3  =  7,  COL4  =42 

CLS 


'  Display  key  instructions 


LOCATE  1,  COL1 

PRINT  "UP  .  Move  to  next  field" 

LOCATE  2,  COL1 

PRINT  "DOWN  .  Move  to  previous  field" 

LOCATE  3,  COL1 

PRINT  "LEFT/RIGHT  ....  Change  field  up/down" 

LOCATE  4,  COL1 

PRINT  USING  "  F2  .  Turn  sound  to  &  flag5 

LOCATE  5,  COL1 

PRINT  "ENTER  ....  Start  with  current  values" 

LOCATE  6,  COL1 

PRINT  "ESCAPE  . . .  Return  to  previous  screen" 


'  Display  fields 

LOCATE  ROW,  COL3 :  PRINT  "Type  of  inputs  (stress  or  strain)"; 
LOCATE  ROW,  COL4 :  PRINT  USING  "[  &  ] " ;  menu2 . inputs ; 


LOCATE 

ROW 

+ 

2, 

COL3  : 

PRINT 

"Equation" ; 

LOCATE 

ROW 

+ 

2, 

C0L4  : 

PRINT 

USING  " t  &  ]  " 

;  equationname 

LOCATE 

ROW 

+ 

4, 

COL3  : 

PRINT 

"Fixed  /  Varing  n'  and  K' " ; 

LOCATE 

ROW 

+ 

4, 

COL4  : 

PRINT 

USING  " [  &  ]  " 

;  nKtype ; 

LOCATE 

ROW 

+ 

6, 

COL  3  : 

PRINT 

"Calculation 

Type  "  ; 

LOCATE 

ROW 

+ 

6, 

COL4  : 

PRINT 

USING  " [  &  ] " 

;  calctype; 

LOCATE 

ROW 

+ 

8, 

COL3  : 

PRINT 

"Input  file  name  "; 

LOCATE 

ROW 

+ 

8, 

COL4  : 

PRINT 

USING  " [  &  ]  " 

;  inputfile; 

LOCATE 

ROW 

+ 

10, 

,  COL3 

:  PRINT  "Output  file 

name  "  ; 

LOCATE 

ROW 

+ 

10, 

,  COL4 

:  PRINT  USING  " [  &  ] 

";  outputfile; 

END  SUB 
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AB.  FUNCTION  xxKf 


===========================  xxKf  ================================= 

Returns  the  value  calculated  for  K'  based  on  the 

number  of  cycles  executed  < 


FUNCTION  xxKf 

SHARED  menu2  AS  Options2 


SELECT  CASE  menu2 . option3 
CASE  1  ' fixed  K' 

xxKf  =  Kf 

CASE  2  '  variable  K' 

IF  NNf count  >  1000000  THEN 
xxKf  =  Kf 

ELSEIF  NNf count  <  2000  THEN 
xxKf  =  K 
ELSE 

xxKf  =  ( (Kf  -  K)  /  (LOG10  (500000)  -  LOG10 (1000) ) )  * 

(LOG10 (NNf count  /  2)  -  LOGIO(IOOO))  +  K 

END  IF 

CASE  3  'fixed  K' 

xxKf  =  Kf 

CASE  4  '  variable  K' 

IF  NNf count  >  1000000  THEN 
xxKf  =  Kf 

ELSEIF  NNf count  <  2000  THEN 
xxKf  =  K 
ELSE 

xxKf  =  ((Kf  -  K)  /  (LOG10 (500000)  -  LOG10 (1000) ) )  * 

(LOG10 (NNf count  /  2)  -  LOGIO(IOOO))  +  K 

END  IF 
CASE  ELSE 
END  SELECT 

END  FUNCTION 


AC.  FUNCTION  xxnf 


'  =============================  xxnf  ============= 

'  Returns  the  value  calculated  for  n'  based  on  the 
'  number  of  cycles  executed 

FUNCTION  xxnf 

SHARED  menu2  AS  Opt ions 2 

SELECT  CASE  menu2 . option3 

CASE  1  '  fixed  n' 

xxnf  =  nf 

CASE  2  '  variable  n' 

IF  NNf count  >  1000000  THEN 
xxnf  =  nf 
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ELSEIF  NNf count  <  2000  THEN 
xxnf  =  n 
ELSE 

xxnf  =  ( (nf  -  n)  /  (LOGIO ( 500000 )  -  LOGIO (1000) ) )  * 

(LOG10 (NNf count  /  2)  -  LOGIO (1000))  +  n 

END  IF 

CASE  3  '  variable  n' 

IF  NNf count  >  1000000  THEN 
xxnf  =  nf 

ELSEIF  NNf count  <  2000  THEN 
xxnf  =  n 
ELSE 

xxnf  =  ( (nf  -  n)  /  (LOGIO (500000)  -  LOGIO (1000) ) )  * 

(LOGIO (NNf count  /  2)  -  LOGIO (1000))  +  n 

END  IF 

CASS  4  '  fixed  n' 

xxnf  =  nf 
CASE  ELSE 
END  SELECT 

END  FUNCTION 


AD.  FUNCTION  xxNNf count 


'  Returns  the  number  of  cycles  executed 

FUNCTION  xxNNf count & 

'DIM  xxNNf count  AS  LONG 

xxNNf count  =  (top  *  block  +  i) 

END  FUNCTION 
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APPENDIX  B.  PROGRAM  LOADGEN 


Random  load  generation  program 

based  on  a  typicical  1000  hour  block  for  an  A- 6  aircraft 


DIM  outfile  AS  STRING 
CONST  fourg  =  1978 
CONST  fiveg  =  333 
CONST  sixg  =48 
CONST  seveng  =  10 

CONST  totalg  =  fourg  +  fiveg  +  sixg  +  seveng 
CONST  Su  =  84 
CONST  Sy  =  76 
CONST  gdesign  =  6.5 


RANDOMIZE  TIMER 
CLS 

'LOCATE  20,  10:  PRINT  "Enter  name  for  stress  load  output  file:" 

'LOCATE  20,  50:  INPUT  outfile 

CLS 

FOR  j  =  l  to  4  '  loop  to  create  four  random  files 

four count  =  fourg 
five count  =  fiveg 
sixcount  =  sixg 
sevencount  =  seveng 

totalcount  =  fourg  +  fiveg  +  sixg  +  seveng 

SELECT  CASE  j 
CASE  1 


OPEN 

"testaa" 

FOR 

OUTPUT 

AS 

#2 

CASE  2 
OPEN 

"testbb" 

FOR 

OUTPUT 

AS 

#2 

CASE  3 
OPEN 

"testcc" 

FOR 

OUTPUT 

AS 

#2 

CASE  4 
OPEN 

"testdd" 

FOR 

OUTPUT 

AS 

#2 

CASE  ELSE 


END  SELECT 
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g"  load  history  greneration: 


OPEN  "gseries "  FOR  OUTPUT  AS  #1 

WHILE  totalcount  >  0 
DO 

x  =  RND 
xx  =  RND 

IF  x  <=  (seveng  /  totalg)  THEN 
'  7+  "g"  case 

IF  sevencount  =  0  THEN  EXIT  DO 
y  =  7  +  ( INT (xx  *  10) )  /  10 
WRITE  #1,  y 
WRITE  #1,  1 

sevencount  =  sevencount  -  l 
ELSEIF  x  <=  ( (seveng  +  sixg)  /  totalg)  THEN 
'  6  to  7  "g"  case 
IF  sixcount  =  0  THEN  EXIT  DO 
y  =  6  +  ( INT (xx  *  10) )  /  10 
WRITE  #1,  y 
WRITE  #1,  1 

sixcount  =  sixcount  -  1 

ELSEIF  x  <=  ( (seveng  +  sixg  +  fiveg)  /  totalg)  THEN 
'  5  to  6  "g"  case 
IF  fivecount  =  0  THEN  EXIT  DO 
y  =  5  +  ( INT (xx  *  10) )  /  10 
WRITE  #1,  y 
WRITE  #1,  1 

fivecount  =  fivecount  -  1 
ELSE 

'  4  to  5  "g"  case 

IF  four count  =  0  THEN  EXIT  DO 

y  =  4  +  (INT  (xx  *  10)  )  /  10 

WRITE  #1,  y 

WRITE  #1,  1 

fourcount  =  fourcount  -  l 
END  IF 


totalcount  =  fourcount  +  fivecount  +  sixcount  +  sevencount 

LOOP 

WEND 

WRITE  #1,  999 
CLOSE  #1 
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'  conversion  of  "g"  load  history  to  a  stress  load  history: 

IF  Su  >  (1.5  *  Sy)  THEN 

gtostress  =  Su  /  (1.5  *  gdesign) 

ELSE 

gtostress  =  Sy  /  (gdesign) 

END  IF 

OPEN  "gseries"  FOR  INPUT  AS  #3 

INPUT  #3,  load 

WHILE  load  <  999 

stressload  =  load  *  gtostress 
WRITE  #2,  stressload 
INPUT  #3,  load 
WEND 

CLOSE  #2 
CLOSE  #3 

load  =  0 

NEXT  j 

END 
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APPENDIX  C.  MATLAB  DATA  REDUCTION  PROGRAM 


A.  PROGRAM  CODE 

%  UNCORRECTED  STRAIN  +  LOAD  DATA 

format  compact 


%  Monotonic  #1 
%  (81  points) 


%  0 

1 

2 

3 

4 

exOa=  [ . 0004 

.0012 

.  0022 

.  0030 

.0038 ; 

.0091 

.0099 

.0110 

.0119 

.0127; 

.  0176 

.  0184 

.  0192 

.  0200 

.0208; 

.  0257 

.0265 

.  0273 

.0281 

.0289 ; 

.  0356 

.  0381 

.  0401 

.  0422 

.  044  3  ; 

.  0573 

.  0595 

.  0617 

.  0640 

.0663  ; 

.  0859 

.  0880 

.0901 

.0923 

.  0944 ; 

.2005 

.2376 

.2727 

.3034 

.3316; 

.5399 

.5786] 

/ 

ld0a=[  305 

628;  998 

1332 

1650; 

3724 

4118 

4488 

4829 

5145  ; 

7044 

7345 

7663 

7973 

8278; 

10140 

10448 

10726 

11036 

11341; 

13738 

14374 

14874 

15312 

15695; 

16971 

17098 

17214 

17307 

17390; 

17884 

17921 

17967 

17990 

18027; 

18901 

18843 

18916 

18992 

19065; 

19326 

19339]  ; 

%  "onotonic  #2 

%  (88  points) 

%  0 

1 

2 

3 

4 

ex0b=  [ . 0005 

.0012 

.0021 

.0029 

.0037; 

.  0085 

.0093 

.0103 

.  0112 

.0120; 

.0166 

.0174 

.0182 

.0190 

.  0198 ; 

.0244 

.  0251 

.0259 

.0267 

.0274; 

.0340 

.0362 

.0382 

.  0404 

.0425; 

.0553 

.  0574 

.0596 

.0618 

.0641; 

.0812 

.0834 

.0855 

.  0877 

.  0900 ; 

.  1934 

.2134 

.2336 

.2538 

.  274  2 ; 

.3958 

.4159 

.4359 

.4561 

.4762  ; 

ld0b=[  306 

619 

962;  1303 

1619  ; 

3596 

3972 

4375 

4697 

5001; 

6890 

7195 

7511 

7815 

8131; 

9934 

10240 

10540;  10852 

11134  ; 

13536 

14209 

14795 

15374 

15854  ; 

17399 

17496 

17586 

17667 

17744 ; 

18095 

18131 

18163;  18197 

18221; 

18851 

18927 

19002 

19052 

19105; 

19224 

19258;  19381;  19404 

19417; 

(for  K  and  n) 


5  S  7  8  9 


0047 

.  0055 

.0063 

.  0071 

.  0081 

0136 

.  0143 

.  0152 

.  0159 

.0168 

0216 

.  0224 

.  0232 

.  0241 

.  0249 

0297 

.0305 

.0313 

.  0322 

.  0330 

0464 

.  0487 

.  0510 

.0530 

.  0551 

0687 

.  0707 

.  0721 

.0756 

.0835 

0966 

.  1170 

.  1373 

.  1574 

.  1777 

3604 

.3924 

.4274 

.4643 

.5019 

1979; 

2307; 

2627;  2975;  3350 

54  63  ; 

5771  ; 

6088;  6406;  6731 

8574; 

8896; 

9214;  9526;  9841 

11635; 

11943; 

12223;  12513;  12789 

16003; 

16278; 

16515;  16687;  16839 

17475; 

17534  ; 

17616;  17675;  17888 

18060; 

18266; 

18419;  18532;  18627 

19128; 

19176; 

19225;  19264;  19302 

5 

6 

7  8  9 

.0045 

.0053 

.0060 

.0068 

.0076 

.  0128 

.0135 

.0143 

.  0151 

.0159 

.  0205 

.  0213 

.  0220 

.0228 

.0236 

.  0282 

.0290 

.0298 

.0305 

.0313 

.0446 

.0468 

.0489 

.0510 

.0531 

.0664 

.0684 

.0704 

.0725 

.0791 

.0923 

.1125 

.1327 

.1530 

.1734 

.2943 

.  3145 

.3350 

.3550 

.3755 

.4965 

.5166 

.5368 

.5572] 

/ 

1958 

2275 

2590 

2911 

3242 

5333 

5631 

5943 

6259 

6569 

8417 

8723 

9041 

9339 

9651 

11442 

11749 

12031 

12300 

12574 

16297 

16676 

16931 

17147 

17273 

17804 

17854 

17901 

17947 

18073 

18248 

18432 

18571 

18661 

18768 

19153 

19200 

19241 

19277 

19305 

19433 

19438 

19449 

19452] 

f 

%  10  percent:  3790  cycles 

%  (73  points) 


%  0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

exl0  =  [ . 0005 

.  0013 

.  0022 

.  0030 

.  0038 

.  0047 

.  0054 

.  0063 

.  0071 

.  0079 

.0088 

.  0099 

.  0108 

.0117 

.  0126 

.  0134 

.  0142 

.  0151 

.  0159 

.  0167 

.  0175 

.  0184 

.  0192 

.0200 

.  0209 

.0217 

.0225 

.  0233 

.  0242 

.  0250 

.  0258 

.0266 

.  0274 

.  0283 

.  0291 

.0299 

.  0307 

.  0316 

.  0324 

.  0332 

.  0361 

.  0381 

.  0402 

.  0422 

.  0444 

.  0466 

.  0489 

.  0511 

.0532 

.  0553 

.  0574 

.  0595 

.0618 

.0638 

.  0660 

.  0682 

.  0703 

.  0725 

.  0746 

.  0778 

.  0798 
.  1911 

.  0818 
.  2114 

.  0840 
.2314 

.  0861 
.2516: 

.  0882 

.0902 

.  1102 

.  1304 

.  1506 

.  1707 

ldl0=[  302 

595 

934 

1252 

1561 

1864 

2160 

2469 

2770 

3098 

3446 

3795 

4170 

4486 

4785 

5076 

5365 

5671 

5967 

6274 

6571 

6856 

7161 

7450 

7748 

8026 

8334 

8619 

8911 

9211 

9490 

9787 

10071 

10361 

10652 

10946 

11225 

11517 

11812 

12074 

13059 

13751 

14478 

15148 

15875 

16591 

17315 

18004 

18577 

19073 

19372 

19536 

19560 

19573 

19577 

19572 

19553 

19524 

19483 

19365 

19325 

19202 

19297 

19137 

19303 

19015 

19297 

18802; 

19282 

/ 

19324 

19304 

19264 

19261 

19243 

%  20  percent:  7580  cycles 

%  (87  points) 


%  0 
°r 

1 

2 

3 

4 

5 

6 

7 

8 

9 

■q - —  —  —  — •  —  —  —  — 

ex20= [ . 0007 

.0016 

.  0028 

.0039 

.  0049 

.  0060 

.  0069 

.  0081 

.  0090 

.0102 

.  0113 

.  0121 

.0132 

.  0142 

.0153 

.0162 

.  0172 

.0181 

.0191 

.  0201 

.  0211 

.0220 

.  0229 

.0240 

.0250 

.0260 

.0271 

.0282 

.0295 

.0309 

.0322 

.0334 

.0349 

.  0363 

.0379 

.0394 

.  0410 

.  0426 

.  0442 

.0458 

.0515 

.  0536 

.0558 

.0581 

.0604 

.  0626 

.0649 

.  0672 

.  0695 

.  0717 

.  0740 

.  0762 

.  0784 

.  0806 

.  0826 

.0849 

.  0873 

.  0897 

.  0918 

.  1010 

.  1031 

.  1051 

.  1073 

.  1094 

.  1117 

.  1139 

.  1339 

.  1542 

.  1747 

.  1948 

.2150 

.2351 

.  2553 

.2756 

.  2959 

.3159 

.3360 

.3561 

.3764 

.  3965 

.4168 

.4368 

.4569 

.4771 

.4974 

.  5179 

.5385 

.5586] 

• 

ld20= [  327 

650 

984 

1324 

1641 

1955 

2271 

2609 

2952 

3291 

3592 

3911 

4246 

4540 

4896 

5206 

5513 

5827 

6132 

6468 

6782 

7077 

7411 

7699 

8034 

8301 

8629 

8918 

9245 

9516 

9816 

10093 

10364 

10670 

10939 

11241 

11518 

11793 

12051 

12331 

13238 

13543 

13876 

14198 

14513 

14846 

15122 

15467 

15741 

16030 

16302 

16577 

16838 

17029 

17194 

i  7366 

17483 

17593 

17656 

17921 

17944 

18007 

18062 

18086 

18120 

18157 

18402 

18583 

18703 

18802 

18888 

18951 

19025 

19092 

19140 

19209 

19253 

19297 

19343 

19387 

19414 

19456 

19463 

19501 

19528 

19576 

19567 

19571] 

; 

%  30  percent:  11370  cycles 

%  (88  points) 

%  0  12  3 

4 

5 

6 

7 

8 

9 

ex30= [ . 0004 

.  0012 

.  0021 

.  0030 

.  0038 

.  0046 

.  0055 

.0063 

.  0071 

.  0079 

.  0089 

.0098 

.  0107 

.  0116 

.  0125 

.0133 

.  0141 

.  0150 

.  0158 

.  0166 

.  0174 

.  0182 

.  0191 

.  0199 

.  0207 

.  0215 

.  0223 

.  0230 

.  0239 

.  0246 

.0255 

.  0263 

.0271 

.0279 

.  0278 

.0295 

.  0303 

.  0311 

.  0319 

.  0327 

.  0355 

.  0376 

.0398 

.0419 

.  0441 

.  0463 

.  0483 

.  0505 

.  0527 

.  0547 

.  0568 

.  0589 

.0609 

.0630 

.0651 

.0693 

.  0716 

.  0736 

.  0819 

.  0843 

.0843 

.  0867 

.0892 

.  0915 

.0938 

.0963 

.  1169 

.1372 

.  1576 

.  1780 

.1989 

.2196 

.2397 

.2602 

.2803 

.3008 

.3209 

.3412 

.3613 

.3813 

.4016 

.4216 

.4418 

.4625 

.4829 

.5031 

.5235 

.5438 

.5645] 

' 

ld3  0=  [  293 

595 

939 

1266 

1583 

1912 

2225 

2539 

2871 

3196 

3546 

3901 

4302 

4632 

4917 

5237 

5562 

5871 

6176 

6486 

6812 

7112 

7423 

7719 

8029 

8321 

8608 

8900 

9224 

9520 

9809 

10117 

10402 

10719 

11009 

11293 

11599 

11905 

12193 

12473 

13473 

14261 

15044 

15811 

16554 

17289 

17942 

18590 

18954 

19085 

19101 

19087 

19093 

19094 

19091 

19085 

19086 

19105 

19098 

19089 

19093 

19073 

19114 

19104 

19111 

19088 

19098 

19053 

19081 

19132 

19181 

19239 

19255 

19295 

19297 

19332 

19365 

19389 

19438 

19486 

19489 

19498 

19501 

19493 

19498 

19487 

19503 

19510 

19498] 

; 

%  40  percent:  15160  cycles 


%  (86  J 

%  0 

joints) 

1 

2 

3 

4 

5 

6 

7 

8 

9 

0 - 

ex40= [ .0002 

.0015 

.0031 

.0045 

.0057 

.0070 

.0082 

.0095 

.0107 

.0119 

.0129 

.0144 

.0158 

.0173 

.0187 

.0205 

.0220 

.0236 

.0253 

.0269 

.0286 

.0304 

.0320 

.0337 

.0354 

.0372 

.  0389 

.  0408 

.0425 

.  0442 

.  0460 

.0479 

.  0497 

.0515 

.0534 

.  0554 

.0575 

.0593 

.  0612 

.  0631 

.0697 

.0717 

.  0738 

.  0759 

.0780 

.  0800 

.0820 

.  0845 

.0866 

.  0891 

.  0911 

.0934 

.0958 

.0983 

.  1005 

.1028 

.1051 

.1075 

.  1097 

.  1191 

.  1214 

.  1237 

.  1261 

.1283 

.  1306 

.1328 

.1532 

.  1737 

.  1939 

.2139 

.2343 

.2545 

.  2746 

.  2950 

.3151 

.3355 

.3555 

.  3758 

.3962 

.4165 

.4369 

.4571 

.4772 

.4973 

.  5175 

.5379 

.5581] 

• 

ld40=[  158 

450 

835 

948 

1230 

1682 

2094 

2415 

2720 

2949 

3194 

3482 

3794 

4049 

4339 

4635 

4942 

5220 

5500 

5777 

6077 

6356 

6625 

6915 

7192 

7509 

7741 

8026 

8301 

8564 

8836 

9103 

9381 

9644 

7923 

9507 

9712 

10160 

10471 

10657 

9823 

12246 

10089 

10914 

12698 

11817 

10769 

13885 

13498 

11450 

12099 

13279 

14430 

14394 

14704 

15653 

14991 

16104 

16423 

16972 

17423 

16899 

17247 

17428 

17435 

17558 

18103 

18422 

1861' 

18789 

18920 

19008 

19074 

19148 

19206 

19185 

19232 

19362 

19395 

19464 

19503 

19539 

19445 

19588 

19634 

19664 

19667] 
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%  adjust  number: 


area  =  .046875; 


ldOa  =  ldOa 

/ 

5 

stressOa  =  ldOa 

/ 

area; 

ldOb  =  ldOb 

/ 

5 

stressOb  =  ldOb 

/ 

area  ; 

ldlO  =  ldlO 

/ 

5 

stresslO  =  ldlO 

/ 

area; 

ld20  =  ld20 

/ 

5 

stress20  =  ld20 

/ 

area; 

ld30  -  ld30 

/ 

5 

stress30  =  ld30 

/ 

area; 

ld40  =  ld40 

/ 

5 

stress40  =  ld40 

/ 

area ; 

strainOa 

= 

exOa 

★ 

.  3 

strainOb 

s 

exOb 

* 

.  3 

strainlO 

= 

exlO 

* 

.3 

strain20 

= 

ex20 

★ 

.3 

strain30 

ex30 

★ 

.3 

strain40 

= 

ex40 

★ 

.  3 

%!del  a:\ernie.met 
%!del  ermie.met 


plot (strainOa, stressOa) , title ( ' Monotonic  #1  Stress-Strain' ) , grid 
xlabel (' Strain  (in/ in) ') ,ylabel (' Stress  (psi)'); 

%pause 

%meta  a:\ernie 


plot (strainOb, stressOb) , title ( ' Monotonic  #2  Stress -Strain' ) , grid 
xlabel (' Strain  (in/in) ') .ylabel (' Stress  (psi)'); 

%pause 

%meta  a:\emie 

plot (strainlO , stresslO) , title (' 10%  Life  Stress-Strain' ), grid 
xlabel (' Strain  ( in/in) '), ylabel (' Stress  (psi)'); 

%pause 

%meta  a:\ernie 

plot (strain20 , stress20) , title (' 20%  Life  Stress-Strain' ), grid 
xlabel (' Strain  (in/in) '), ylabel (’ Stress  (psi)'); 

%pause 

%meta  a:\ernie 


plot (strain30, stress30) , title (' 30%  Life  Stress-Strain' ), grid 
xlabel (' Strain  (in/in) ') .ylabel (' Stress  (psi)'); 

%pause 

%meta  a:\ernie 

plot (strain40 , stress40) , title (' 40%  Life  Stress-Strain' ) .grid 
xlabel (' Strain  (in/in) ') .ylabel ('Stress  (psi)'); 

%pause 

%meta  a:\ernie 


plot (strainOb, stressOb, strainlO , stresslO , strain20, stress20 , . . . 
strain30, stress30, strain40, stress40) , title ( ' Stress -Strain' ) , grid 
xlabel (' Strain  (in/in) ') .ylabel (' Stress  (psi)'); 

%pause 

%meta  a:\emie 


76 


!del  ernie.out 
diary  erine.out 


i 


%  compute  true  strain  from  engineering 


tstrainOa  =  log 
tstrainOb  =  log 
tstrainlO  =  log 
tstrain20  =  log 
tstrain30  =  log 
tstrain40  =  log 


(1 

+ 

strainOa) 

(1 

+ 

strainOb) 

(1 

+ 

strainlO) 

(1 

+ 

strain20) 

(1 

+ 

strain30) 

(1 

+ 

strain40) 

%  compute  true  stress  from  engineering 


sigl  =  stressOa  . * 
sig2  =  stressOb  . * 
sig3  =  stresslO  .* 
sig4  =  stress20  .* 
sig5  =  stress30  .  * 
sig6  =  stress40  .* 


(l  +  tstrainOa) 
(1  +  tstrainOb) 
(1  +  tstrainlO) 
(1  +  tstrain20) 
(l  +  tstrain30) 
(1  +  tstrain40) 


lsigl  =  loglO (sigl); 
lsig2  =  loglO (sig2)  ; 
lsig3  =  loglO (sig3); 
lsig4  =  loglO (sig4) ; 
lsig5  =  loglO (sig5); 
lsig6  =  loglO (sig6) ; 


E  =  10300000; 


xl  =  loglO (tstrainOa (2 : 82)  -  sigl (2 : 82) /E) 

x2  =  loglO (tstrainOb  -  sig2/E) ; 

x3  =  loglO (tstrainlO  -  sig3/E) ; 

x4  =  loglO (tstrain20  -  sig4/E) ; 

x5  =  loglO (tstrain30 (2 : 89)  -  sig5 (2 : 89) /E) 

x6  =  loglO (tstrain40 (2 : 87)  -  sig6 (2 : 87) /E) 


lsigla  =  lsigl (2 : 82) ; 
lsig5a  =  lsig5 (2:89)  ; 
lsig6a  =  lsig6 (2 : 87) ; 


xxl 

xx2 

xx3 

xx4 

xx5 

xx6 


polyf it (xl , 
polyf it (x2 , 
polyf it (x3 , 
polyf it (x4 , 
polyf it (x5 , 
polyf it (x6 , 


lsigla,  1) 
lsig2,  1) 
lsig3 ,  1) 
lsig4,  1) 
lsig5a,  1) 
lsig6a,  1) 
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K1 

= 

10* (xxl (2) ) 

K2 

= 

10* (xx2 (2) ) 

K3 

= 

10* (xx3 (2) ) 

K4 

= 

10* (xx4 (2) ) 

K5 

s 

10* (xxS (2) ) 

K6 

= 

10* (xx6  (2) ) 

nl 

= 

xxl (1) 

n2 

= 

xx2 ( 1 ) 

n3 

= 

xx3  (1) 

n4 

= 

xx4 (1) 

n5 

= 

xx5 (1) 

n6 

= 

xx 6 (1) 

diary 

B.  PROGRAM  OUTPUT/RESULTS 


xxl  = 

0.4425 

5.6106 

xx2  = 

0.4331 

5.5810 

xx3  = 

0.6246 

6.0503 

xx4  = 

0.5063 

5.6425 

xx  5  = 

0.4261 

5.5662 

xx  6  = 

0.5902 

5.6822 

K1  = 

4 . 0798e+005 

K2  = 

3 . 8107e+005 

K3  = 

1. 1229e  +  006 

K4  = 

4 . 3900e  +  005 

K5  = 

3 . 6831e+005 

K6  = 

4 . 8105e+005 

nl  = 

0.4425 

n2  = 

0.4331 

n3  = 

0.6246 

n4  = 

0.5063 

n5  = 

0.4261 

n6  = 

0.5902 

APPENDIX  D.  MATERIAL  DATA  BASE 


A.  BRITISH/AMERICAN  UNITS 


Monotonic  and  Cyclic  Strain  Properties  ol  Selected  Engineering  Alloys:  American/British  unite 


Material 

Process 

Descnpiion 

.s*. 

(ksit 

S  |f  S  V 
tksirksii 

KlK‘ 

\  ksi/k^i) 

nm ' 

(fit, 

<r  t  a , 

tksi/ksii 

h 

f 

5, 

(2.V  =7  10’). 

tksi) 

.5,5, 

Steel 

1005-  1009 

H  R  sheet 

50 

IR'17 

77  67 

I)  16  0  12 

1.60.10 

1 2?'q^ 

-0.109 

-O  ?Q 

:i 

0  4' 

11X15- 1009 

C  D  sheet 

60 

58  16 

7671 

0.0490  I) 

1.02/0.11 

12278 

-0.071 

-0.41 

28 

0  47 

1020 

H  R.  sheet 

64 

1815 

107  112 

0  19  0  18 

0.96  0  4 1 

103/130 

-0.12 

-0.51 

22 

0  14 

00.10' 

Casl  steel 

72 

44/4(1 

0.30  0  |i 

0  62/0.28 

109/95 

-0.082 

-0.51 

28 

0.38 

Man- 1  en 

HR  sheet 

74 

57/54 

—  114 

0.20/0.11 

1  02/0  86 

118/117 

-0.071 

-0.65 

?8 

0  51 

1040 

As  forged 

90 

50'56 

— 

0.22/0.18 

0.9.1/0.61 

152/22? 

-0.14 

-0.57 

2< 

0.28 

RQC  100 

H  R  sheel 

13' 

128'87 

170208 

0.06/0.14 

1  02/0.66 

193/580 

-0  07 

-0  69 

59 

0  41 

4142 

Drawn  at  icmp 

154 

152/108 

- - - 

—  0.18 

0.35/0.22 

162/210 

-0  10 

-0  51 

44 

0.28 

4142 

Q  &  T 

20' 

200  120 

— .  — 

0  051/0.17 

0. 66/(1  45 

265  26' 

-0.08 

-0.75 

7? 

0  ?6 

4142 

0  &  7 

280 

250  19' 

— 

0.048/0  13 

0.43/0.09 

315  31' 

-0  081 

-0  61 

8' 

0  1! 

4. 140 

II.R.  and  annealed 

120 

92 '66 

— 

—  '0  18 

0.57/0  45 

158/174 

-0  095 

-0.54 

40 

0  I! 

414(1 

V  &  T 

ISO 

no  no 

0  066/0  14 

0  84/0.73 

240740 

-0.076 

-0  62 

7| 

0  40 

4140 

2H 

|9Q  pH 

— 

—  0  15 

0  48/(1  48 

226790 

-0091 

-0  60 

68 

0 1: 

9262 

Annealed 

H4 

66  76 

251  200 

0.22/0  15 

0  16/0  16 

151151 

-0  071 

-0  47 

50 

0  ?S 

9262 

<)  <ft  I 

145 

114 '94 

—  19' 

0  140  12 

0  41,041 

177/17’ 

-0.071 

-0  60 

55 

0  is 

Aluminum 

1 100-0 

As  received 

Ift 

149 

— 

—'0  1 5 

2  09.1  8 

—  28 

-0  106 

-0  6q 

5 

0  71 

2024-Ti 

— 

68 

"'62 

66'  05 

0.032/0  06' 

0.28 '0.22 

81/160 

-0.124 

-0  59 

0  32 

2024-74 

— 

69 

4J7J 

117— 

0  20  0  08 

0  41.0.21 

92147 

-0  II 

-it?: 

25 

0  r 

5456- H? 

— 

58 

14'52 

—  — 

—70.16 

0  42/0.46 

76/105 

-0  II 

-0  67 

18 

o  11 

7075-7  ft 

— 

84 

68,76 

120- 

0.1 1/0  146 

0.41/0  19 

I0R/ 19 1 

-0  126 

-0.52 

2s 

0  ?0 
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B 


SI  UNITS 


Monotonic  and  Cyclic  Strain  Properties  of  Selected  Engineering  Alloys:  SI  Units 


Material 

Process 

Description 

5.  W- 

iMPal  iMP.xMPai 

K  K' 

i  MPa.- MPa  l 

nn ' 

<Tf  r r,  . 

i  M  Pa.  M  fa ) 

h 

2.V  -=  !0Tt 
•  M  Pa  i 

5 , 5. 

WW 

1005- 1009 

HR.  sheet 

us 

202.228 

5  4 1 4/>: 

0  16(1  12 

!  Or()  10 

848/04  1 

M  |09 

-0  19 

148 

0  41 

1005- 1009 

C.D.  sheet 

1 14 

400/248 

524.290 

n.tu9/t).  1 1 

1  02/0  1 1 

84  1  538 

-0  073 

-0  41 

195 

0.47 

1020 

HR  sheet 

441 

262  241 

U  8.772 

o  190  |8 

0  96,0  41 

-10, 896 

0  12 

-0  51 

152 

0  14 

oo.W 

Cast  steel 

496 

>03  317 

— 

0  30/0  1 1 

0  62.  0  28 

“SO/653 

-0  082 

•0  <1 

190 

0  38 

Man-  Ten 

HR.  sheet 

MO 

191/372 

_780 

0  20/0  1 1 

1  02  0  86 

8|4  807 

o  071 

■  o  o5 

262 

0  51 

1040 

As  forged 

621 

US '380 

— 

0.22/0  18 

0  93/0  01 

1050  1540 

-«i  14 

-it  57 

173 

o  28 

RQC-100 

H  R.  sheet 

93 1 

883/6(8) 

1  172  1434 

0.00/0. 14 

1.02.  0  66 

1  330  1240 

-0  07 

-n  09 

401 

0  41 

4142 

Drawn  at  temp 

1062 

1048/74  s 

— 

—0.18 

0  35  0  22 

1115/1450 

-0  10 

-0  51 

310 

o  28 

4142 

0  A  I 

1413 

1379  827 

—  — 

0.051/0  17 

1"  66,0  45 

1825/1825 

-0  08 

-0  7S 

501 

o  30 

4142 

Q  A  r 

1931 

1724,1344 

— 

0.048/0  13 

0  4 3  0  09 

2170/2170 

-0  081 

-o  01 

'89 

o  31 

4540 

H  R.  and  annealed 

S27 

634.455 

— 

—  0  18 

0  57/0.45 

1090/1200 

-0.095 

-0  S4 

274 

n  33 

4340 

<)  &  r 

1211 

II727S8 

|S7<*_ 

04)66/0  14 

0  84,0  73 

1655/1655 

-i)  076 

o  02 

492 

0  40 

4340 

0  A  T 

1469 

1 472/827 

— 

—0  15 

0  48/0  48 

1560.2000 

-0  091 

o  00 

407 

o  32 

9262 

Annealed 

924 

4SS  524 

1744.1  379 

0  22  0  15 

0  100  10 

1040- 1040 

-o  07 1 

-0  47 

348 

n  38 

9262 

<)  A  r 

1000 

"86648 

—  1358 

0  14/0.12 

0  4  I  /O  4| 

1220  1220 

-9  073 

-o  oo 

>S  1 

n  18 

Mttnunum 

1 100-0 

As  received 

110 

97,62 

— 

—  0  15 

2  09' 1.8 

—  193 

-0  100 

0  09 

37 

it  33 

2024-T3 

— 

469 

>79/427 

455  655 

0.032/0  065 

0  28  0  22 

558/1  1(8) 

-0  124 

-n  59 

151 

o  32 

2024-T4 

— 

476 

'03  441 

807’— 

0.20  0  08 

0  43  o  ;i 

634/1015 

-0  II 

-o  <2 

175 

it  17 

5456-  H. 1 

— 

400 

234/159 

— 

—  0  16 

o  42/0  46 

*24'  725 

oil 

o  07 

124 

0  31 

7075-T6 

— 

579 

469/524 

<?2?' _ 

0  111)  140 

t>  ,41'tV  19 

'45/1315 

-0  120 

o  52 

|-6 

It  10 
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